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Energetic galaxy-wide outflows in high-redshift ultra-luminous 
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ABSTRACT 

We present integral field spectroscopy observations, covering the [O lll]A,A4959,5007 
emission-line doublet of eight high-redshift (z= 1.4-3. 4) ultra-luminous infrared galaxies 
(ULIRGs) that host Active Galactic Nuclei (AGN) activity, including known sub-millimetre 
■ luminous galaxies (SMGs). The targets have moderate radio luminosities that are typical of 

high-redshift ULIRGs (Li. 4 ghz=10 24 -10 25 WHz *) and therefore are not radio-loud AGN. 
We de-couple kinematic components due to the galaxy dynamics and mergers from those due 
to outflows. We find evidence in the four most luminous systems (Lpm] ^ 10 43 ergs -1 ) for the 
signatures of large-scale energetic outflows: extremely broad [O III] emission (FWHM«700- 
1400 kms~') across «4-15kpc, with high velocity offsets from the systemic redshifts (up 
to «850kms _1 ). The four less luminous systems have lower quality data displaying weaker 
evidence for spatially extended outflows. We estimate that these outflows are potentially de- 
positing energy into their host galaxies at considerable rates (E «10 43 -10 45 erg s~'); however, 
due to the lack of constraints on the density of the outflowing material and the structure of the 
outflow, these estimates should be taken as illustrative only. Based on the measured maximum 
velocities (v max « 400- 1400 km s~') the outflows observed are likely to unbind some fraction 
of the gas from their host galaxies, but are unlikely to completely remove gas from the galaxy 
haloes. By using a combination of energetic arguments and a comparison to ULIRGs with- 
out clear evidence for AGN activity, we show that the AGN activity could be the dominant 
power source for driving all of the observed outflows, although star formation may also play 
a significant role in some of the sources. 

Key words: galaxies: high-redshift, sub millimetre; — galaxies: evolution; — galaxies: star 
formation rates, AGN; — galaxies: individual 



1 INTRODUCTION 

Many of the most successful models of galaxy formation require 
energetic outflows over galaxy scales (i.e. ~l-10kpc) to repro- 
duce properties of local massive galaxies (e.g., black-hole-spheroid 



mas s relationship; the bright-end of the galaxy luminosity function ; 



. JSilk & R ees 1998; Springel et al. 2005: Pi Matteo et alj200l : 
Hopkins et al. 2006, 2008; Debuhr et all2012h . Observations show 
that galaxy-wide outflows can be powered b y star formation and/or 
active galactic nuclei (AGN) activity (e.g., lHeckmanet"ai]|l99d : 
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ICrenshaw etaUl2003l; iRupke et alj|2005allbl) ; however, most mod- 
els predict that AGN-driven outflows are required in the most mas- 
sive galaxies. These simulated outflows inhibit further black-hole 
accretion and star formation by injecting considerable kinetic en- 
ergy into the interstellar medium (ISM) of the host galaxies. While 
theoretically attractive, these models can only be judged to be suc- 
cessful if supported by observational data. 

AGN-driven outflows are initially launched from the accre- 
tion disk or dusty torus surrounding the black hole (BH), either 
in the form of a radio jet/lobe or a radiatively-driven wind. Radio 
jets/lobes are thought to be most effective at h eating gas in radio- 
loud systems thro ugh the so-called radio mode dCroton et al . 2006; 
iBower et alj|2006l ) and evidence of this is observed in low-redshift 
massive galaxies, where me chanical heating appears to be prevent- 
mg gas from cooling (e.g. iBest et"all l2005, 2006; Smolcic et al.l 



2009; Dani elson et aTll2012h . In contrast, many models implement 
the more rapid quasar mode, where the radiation field produced 
by AGN launches a high-velocity wind which swee ps up gas in 
the ISM and result s in a kpc-scale "outflow" (e.g ., Springel et al. 
20051; iKind l2005l; |Pi Matteo etall 120051 ; I Hopkins et all 120061 ; 



Debuhr et al.l 20121) . These radiatively-driven outflows are pre- 
dicted to be most prevalent at high redshi ft, (z ~ 2), in luminous 
AGN with high accretion rates (e.g., |Pi Matteo et al. 2005|; I Kind 
l2005h . 

Observationally, AGN-driven winds are known to be preva- 
lent in luminous AGN at low and high redshift: X-ray and UV 
absorption-lin e spectroscopy reveal high-velocity winds of up to 
*0. 1c (e.gjReeyes et alj|2003l;lBlustin et alj|2003l:lTrump et al] 



l2006l;lGibson et alj2009l ; lGofford et al.l201 ll;|Page et al.l201 ll) and 
may be a ubiquitous property of all AGN dGangulv & Brothertonl 
2008). However, these high-velocity winds are most likely to 
be produced close to the accretion di sk (< 1 pc scale; e.g., 
Crensh aw et al .l2003l ; lTombesi et al.l2012h . For these winds to have 
a global effect on the host galaxy, they must couple to the ISM and 
drive the gas over galaxy-wide (~ 1-10 kpc) scales. To test the 
impact of AGN-driven outflows on the formation and evolution of 
galaxies, it is therefore necessary to search for kpc-scale energetic 
outflows. 

Spatially resolved spectroscopy provides a particularly direct 
method to search for and characterise kpc-scale outflows. Such ob- 
servations of low-redshift galaxies have identified kpc-scale out- 



flows in both ionised (e . g.. iHolt et al]|2008l; 



Rupke & V eilleuxll2011 



gas (e.g., Feruglio et al] 



Westmoquette et al 



20 id ; lAlatalo et al 



Fu & Stocktor]|2009l ; 
20121) and molecular 



201 ll) . At high red- 



shift, where accretion-related outflows are predicted to be most 
prevalent, integral-field unit (IFU) o bservations of high-redshift 
radio galaxies (HzRGs, at z ~ 2-3; iNesvadba eta"i]|2006l l2008h 
have showed that kpc-scale energetic outflows are present in at 
least a fraction of the high-redshift galaxy population. AGN-driven 
outflows were revealed in these systems through the presence of 
broad (FWHMMOOOkms- 1 ), high- velocity (>300kms" 1 ) and 
spatially extended [O III] emission^ With implied velocities of 
greater than a few hundred kms -1 these outflows could poten- 
tially drive gas out of the host galaxy, inhibiting future BH growth 
and star formation. Recent IFU observations of the [O III], Ha and 



' Emission from [O III] cannot be produced in high density environ- 
ments, such as the broa d-line region (BLR) of AGN, wi thout being colli- 
sionally de-excited (e.g., lOsterbr ock 1989; Robsonl l 19961) . Luminous (L > 
10 42 erg s" 1 ), broad (FWHM > 1000 km" 1 ), and extended (a; a few kilo- 
parsecs) [O III] emission is therefore not associated with the BLR. 



[C II] emission lin es in high-redshift qua s ars also reveal simila r 
outflow signatures dCano-Dfaz et al.ll2012l ; iMaiolino et alj|2012l) . 
Due to the small number of sources and the biased selection of 
the high-redshift galaxies studied with IFUs so far, the prevalence 
and significance of such galaxy-wide outflows remains poorly con- 
strained. 

Potentially the best place to search for AGN-driven out- 
flows are high-redshift far-IR (FIR) luminous galaxies with Lfir > 
1O 12 L,0 (Ultra-Luminous Infra-Red Galaxies - ULIRGs), such 
as sub-millimetre galaxies (SMGs). These systems are the most 
intensely star-forming galaxies known with star formation rates 
(SFRs) of ^100-1000 M yr~' and are thought to represent a 
rapid star-forming phase tha t every massive galax y goes through 
(e.g., ISwinbank etal] 120061; iTacconi et al] l2008h. They have a 
peak space density at z ~ 2-2.5 I Chapman et al] 120041 . 120051 ; 
Wardl ow et al.l201 lh which roughly coincides with t he peak epoch 
of BH growth (e.g.. [ Cattaneo & Bernardij 120031 ; iHopkins et al] 
120071 ; ISilverman et al.ll2008l) . Using optical, mid-infrared and X- 
ray diagnostics, high-redshift ULIRGs are found to have emis- 
sion dominated by starburst activity with ~30% of SMGs also 
revealing signifi cant AGN activity (e.g., Alexander et al. 2003, 
2005, 2008; lchapman et al.l2005l ; llvison et al.l20ld ; lRaffertv et al] 
2011). Indeed, clustering analysis of SMGs and quasars supports 
the id ea that starburst an d AGN activity occur in the same sys- 
tems dffickox et al]|2012h and that SMGs may evolve into mas- 
sive elliptical galaxies by the present day. Therefore, high-redshift 
ULIRGs that host AGN activity could be in transition from a star- 
formation dominated to an AGN dominated phase. This transition 
phase, which is potentially facilitated by kpc-scale outflows, is an 



evoluti onary period requir e d by many mo dels (e. g.. ISanders et al 
1 19881; ISilk& Reesjjl998l; ISpringel et al] 120051; jpi Matteo et al 
booalHopkins et all2006l.l2008l ; see lAlexander & Hickoxl201lf fo~r 

a review). 

In this paper we present Gemini-North Near-Infrared Inte- 
gral Field Spectrometer (NIFS) and VLT Spectrograph for Inte- 
gral Field Observations in the Near Infrared (SINFONI) IFU ob- 
servations of z ~2 ULIRGs with known optical AGN activity. 
These sources have moderate radio luminosities (Li.4GHz = 10 24 - 
10 25 WHz" 1 ) and can therefore be considered to be radio-quiet 
AGN. In §2 we give details of the IFU observations and data re- 
duction, in §4 we present the analysis and results, in §5 we discuss 
the results and their implications for galaxy evolution models and 
in §6 we give our conclusions. We provide background information 
and detailed results on individual sources in Appendix A. We have 
adopted Hq = 71 kms -1 , £Im = 0.27 and £1/^ = 0.73 throughout; in 
this cosmology, 1" corresponds to 8.5 kpc at z = 2.0. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Target selection 

For this study we selected ten radio-detected z = 1.4-3.4 star- 
forming galaxies that host AGN activity to study their spatially 
resolved dynamics. The sample consists of eight SMGs (labelled 
'SMM' with 850/jm flux densities 5g5o > 4mly) and two sources 
that are undetected at 850,1ml (labelled 'RG'). We will collectively 
refer to the sample as high-redshift ULIRGs throughout. We de- 
tected [O III] emission in eight of the sources (see §4 for de- 
tails). All our detected sources are [O III] luminous (Liq jm= 10 42 - 
10 44 erg s~'), and have high [N Il]/Hcc and [O III]/HP emission-line 
ratios, or an AGN-dominated continuum, indicating the presence of 
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Gemini-North NIFS and VLT SINFONI Observations 
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1 


SMM J0217-0503 


02:17:38.68 


-05:03:39.5 


2.021 


207±16 


5.0 


O 


4.8 


0.47 


4.0 


S 


1 


2 


SMM J0302+0006 


03:02:27.73 


+00:06:53.5 


1.405 


217±9 


2.2 


0,M? 


11.4 


0.42 


3.6 


N 


2 


3 
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03:02:58.94 


+00:10:16.3 


2.239 


55±10 


1.7 


o,u 


9.6 


0.38 


3.2 


N 


1 


4 
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2.315 


160±20 


5.3 
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5 
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60 
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0,U,M 


14.4 
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3.2 


N 


1 


6 


SMMJ1235+6215 


12:35:49.41 


+62:15:36.9 


2.199 


93±6 


2.8 


O.X,M? 


7.8 


0.45 


3.8 


N 


2 


7 


SMM J 1237+6203 
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4.6 
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2.5 


N 


1 


8 
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+40:57:34.5 


2.385 
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8.7 
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2.7 


N 


1 


9 


SMMJ2217+0010 
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+00:10:25.1 


(2.610) 


179±19 


7.9 


0,M 


3.6 


0.40 


3.2 


N 


3 


10 


SMM J22 17+00 17 
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+00:17:02.0 


(2.278) 


114±27 


3.7 


O 


9.6 


0.56 


4.7 


S 
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Table 1. Notes : 

Column (1): Source ID. Column (2): Source name. Columns (3)-(4): The co-ordinates corresponding to th e radio positions from: ILedlow et all (2002); 
Swin bank etaLI 12004); Chapman et al. (2005); Miller et al. 12008); Mo rrison et"ai1 fe01p|) ; lrvison et alj b01 ll) ; Arumugam et al. (in prep); R. J.Ivison (priv. 
comm). Column (5) The redshifts based on the narrow [O III] emission lines using the ga laxy-integrated spec tra p resented in Fig.H] For SMM J2 217+0010 and 
SMM J2217+0017 we do not detect [O III] and we therefore quote the Ha redshift from lTakata et al l fe006j) and|Xlaghband-Zadeh eTaT]j2012l) . respectively. 
Column (6): Radio flux densities taken from the same references as for the radio positions. Column (7): Rest-frame radio luminosities usin g the standard k- 
correction and a spectral index of ct=0.8. Column (8): Ev idence for AGN activity in the sources on the basis of: O: optical emission-line r atios JSwinbank et alj 
2004t Frakata et alj|2006l; lAlaghband-Zadeh et al.ll2012t this work); M: excess in the mid-infrared continuum from infrared spectra ( Valiante*et*aTIl2*0*0*7l; 
Menendez-Delmestre et al. 2009; he evidence for IR-AGN activity in the sources with a "?" is tentative); U: rest-fram e UV spectral signatur es jLedlowetall 
2002 1; ISmail et al JI2003I; IChapman et alj|2004l ; [Chapman et alj|2005t) ; X: X-ray observations ^Alexander et alj|2005l) ; B: Broad-line AGN fcwinbank etai] 
2005; Coppin et al. 2008; see Appendix A for more details on the classifications of individual sources). Column (9): The total on-source exposure times for 
the data presented here. Columns (10)-(1 1): The average seeing for the observations, based on standard star observations. Column (12): The instrument used 
for the observations, either SINFONI (S) or NIFS (N). Column (13) A data quality flag as follows: 1: data which have sufficiently high signal-to-noise ratios 
to study the spatially resolved kinematics; 2. data with lower signal-to-noise rat ios and limited sp atially resolved information; 3. undetected in [O III]^.5007 
in this study. * For SMM J0943+4700 we ob served the radio counter-part H6 (Ledlow et al. 2002) and the 1.4 GHz flux density has been corrected for an 
amplification factor of 1 .2 JCowie et alj2002l) . 



AGN activity (e.g. jKewlev et al1l2006t) . with most of the sources 
having additional observational evidence for AGN activity (see Ta- 
ble[T]for details); see Appendix A for background information on 
the individual sources. We provide details of the target selection 
here. 

Seven of our targets are selected from Tak ata et al. ] d2006l) 
who presented near-IR spectroscopy around the redshifted 
[O III]aX4959,5007 emission-line doublet of 22 high-redshift 
ULIRGs. In order to select sources that might contain AGN- 
driven outflows, we chose seven targets with broad (FWHM > 
700 km s _1 ) and bright (F[ OIII ] >few x 10 _16 ergs _1 ) [O III] emis- 
sion. We al so include in our sample three, z ~ 2-2.5 ULIRGs se- 
lected from Alaghband-Zad eh et al. I l l2012l) which were identified 
as hosting AGN activity on the basis of their rest-frame optical 
emission-line ratios. These targets were initially selected on the 
basis of their Ha luminosities and are therefore not pre-selected 
to have luminous and broad [O III] emission. On e of the sources 
was p resented in a pilot study (SMMJ1237+6203; Alexan der et al.1 
|2010|) and we present the whole sample here. 

Details of our observations and the targets' 1.4 GHz flux den- 
sities, obtained from the literature, are provided in Table 1. The 
1.4 GHz radio luminositie s (Li.4g hz) are calculated using Equation 
(2) from lAlexander et al l 12003), with the 1.4 GHz flux densities 
shown here and an assumed spectral index of a=O.80 We estimate 

2 This value is approximately the mean value for radio-identified SMGs 
with an observ ed range of roughly oc=0.2-1.5, which includes AGN 
hbaretalj|20ld) . Here we define a spectral index a such that flux density, 
S v , and frequency, v, are related by 5 V °= v~ a . 



AGN luminosities for our sources using the SED fitting procedures 
outlined in ij3]and these values are provided in Table|2] 

Although this sample is heterogeneous, it is more represen- 
tative of the overall AGN population than previous IFU studies 
searching for AGN-driven outflows in high-redshift radio galax- 
ies (with L lAGHz Z 10 27 WHz" 1 ; e.g jNesvadba et alj|2008h . For 
example the radio luminosities of our sources (Li.4ghz=10 24 - 
10 25 WHz~') imply a space density of <I> ~10 -5 Mpc~ 3 com- 
pared to <ft <10~ 8 Mpc~ 3 for high-red shift radio galaxies (Fig. [T] 
ISimpson et al]|2012l ; Iwiilott et al.|[l998l) . This is important if we 
are to understand the prevalence of AGN-driven outflows in high- 
redshift galaxies. In Fig. [T] we compare the integrated [O III] prop- 
erties and radio luminosities of our sample with low- and high- 
redshift AGN. 



2.2 Gemini-North NIFS observations 

Seven of the so urces in our sample (those selected from 
iTakata et al.ll2006l) were observed using Gemini-North NIFS (Ta- 
ble 1). The NIFS IFU uses an image slicer to reimage a 3.0 x 3.0 
arcsec 2 field into 29 slices of width 0.103" with a pixel scale of 
0.04" along the slices. The dispersed spectra from the slices are 
reformatted on the detector to provide two-dimensional spectro- 
imaging. Depending on the target redshift, we used the /— , 
H— or K— band grisms to ensure that we cover the rest-frame 
[O III]a,A,4959,5007 emission-line doublet. These grisms have an 
approximate spectral resolutions of A./AA, ~ 6040, 5290 and 5290, 
respectively (FWHM w 2, 3, 4 A, or Av « 50, 57, 57 km s -1 ). For 
our anaylsis we measured the instrumental dispersion more accu- 
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Figure 1. FWHM versus total [O III]^.5007 luminosity (left) and FWHM verses 1.4 GHz luminosity (right) for our eight [O III] detected targets. Shown 
for comparison are low redshift (z<0A) AGN from the Sloan Digital Sky Survey (Mullaney et al. in p rep) and z < 0.5 obscured q uasars with long-slit data 
Greene et al 20111). A lso shown are other high-re dshift sources with IF U data: z~2-3 radio galaxies (Nesvadba et al. 2006, 2008), compact radio galaxies 
NesvadbTetaUl2007l a) , zs±3.5 obscured quasars l lNesvadba et alj|201 lb and a z=2.4 quasar jCano-DiaTetani20*12Tr Our targets have very broad [O III] 



emission, comparable with high-redshift radio galaxies but having 3-4 orders of magnitude lower radio luminosities such that they are ts3 orders of magnitude 
more common at z ~2— 3. When the [O III] profile is decomposed into multiple components, the FWHMs of the broadest components have been plotted. All 
rest-frame radio luminosities are calculated using the standard ^-correction, on the basis of their measured 1.4 GHz flux densities. Radio spectral indices are 
taken to be 01=0.8; however, the horizontal error bars ind icate the range of lum inosities using spectral indic es in the range a=0.2-1.5. The arrows in the right 
panel correspond to space densities of <£> «10~ 5 Mpc~ 3 ^Simpson et al.lE012T) and RilO 
atz «2.5. 



8 5 Mpc 3 ( Willott et al. 1998) from the radio luminosity function 



rately, by measuring the widths of several bright sky-lines close to 
the observed wavelengths of the emission lines for the observations 
of each source. We found that the uncertainties in the instrumental 
dispersion (by measuring the scatter from the widths of several sky- 
lines) to be insignificant (~ 5%) compared to the uncertainties from 
our fitting procedures. We consequently do not include the uncer- 
tainties on the instrumental resolution in our measurements. These 
observations were taken between 2008 March 25 and May 19 and 
2009 August 5 and December 30. All observations were taken in 
<0.5" seeing. 

The observations were performed using repetitions of the stan- 
dard ABBA configuration in which we chopped either 6" to blank 
sky or, if the source was compact (emission-line region <1.5"), we 
chopped within the IFU. Individual exposures were 600s. Final on- 
source exposure times and seeing measurements, calculated using 
the corresponding standard star observations, are shown in Table 1. 

We reduced the data with the Gemini IRAF NIFS pipeline 
which includes sky subtraction, wavelength calibration, and flat 
fielding. Attempts to remove residual OH sky emission lin es were 
made using the sky-subtraction techniques described in iDaviej 
d2007l) . We further corrected for variations in the night sky emission 
by masking the source (when clearly detected) and taking a median 
of each pixel-column and consequently subtracting the (empty sky) 
background separately from each wavelength. 

2.3 VLT SINFONI observations 

Three of the sources in our sample were observed using VLT 
SINFONI (Table 1 ). The full details of these SIN FONI observa- 
tions are given in Alaghband-Zad eh et al. I J2012I) with only ba- 



sic information given here. We used the widest available field-of- 
view, 8.0 x 8.0 arcsec 2 , which is divided into 32 slices of width 
0.25" with a pixel scale 0.125" along the slices . The H+K grating 
was used, which gives an approximate resolution of X/AX ~1500 
(FWHM ~ 1 3 A or Av ~ 200 km s~ 1 ). We made accurate measure- 
ments of the instrumental dispersion for the observations for each 
source in the same manner as outlined for the NIFS observations 
in the previous section. Again we found the uncertainties on the in- 
strumental dispersion to be insignificant compared to the uncertain- 
ties from our fitting procedures. The targets were chopped around 
quadrants of the IFU for sky subtraction. All observations were 
taken in <0.6" seeing. Individual exposures were 600s. Final on- 
source exposure times and seeing measurements, calculated using 
the corresponding standard star observations, are shown in Table 1. 

The data were reduced using the SINFONI ESOREX pipeline, 
which includes sky subtraction, wavelength calibration and flat 
fielding. 



2.4 Flux calibration and stacking 

Each night's data were flux calibrated separately using observa- 
tions of standard stars, at a similar airmass to the target galaxies, 
which were reduced in the same manner as the targets. Since the 
[O III]XX4959,5007 emission-line doublets typically lie in regions 
of > 95% sky transparancy, no corrections were made for telluric 
absorption. To find spatial centroids for the flux-calibrated data 
cubes, white light (wavelength collapsed) images around the red- 
shifted [O III]X5007 emission line (NIFS targets) or Ha emission 
line (SINFONI targets) were produced. The data cubes for which 
we could not measure a centroid were discarded to improve the 
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reliability of the spatially resolved data. Data cubes were only dis- 
carded for SMM J0302+0006 (2 discarded out of the 21 observed), 
RG J0302+0010 (4 discarded out of 20) and SMM J1235+6215 (8 
discarded out of 21). We then spatially aligned and co-added the in- 
dividual data cubes to create the final mosaic, using a median with 
a 3o clipping threshold to remove remaining cosmetic defects and 
cosmic rays. The total on-source exposure times used in the final 
stacks are given in Table 1 . 

In all of the following sections, the effect of instrumental dis- 
persion was corrected for by subtracting it in quadrature from the 
observed line dispersions. 



3 STAR-FORMATION RATES AND AGN LUMINOSITIES 

To assist in the analyses of this work we need to constrain the star- 
formation rates and AGN luminosities of the eight [O IIl]-detected 
sources in a consistent manner. To achieve this we obtained infrared 
flux densities from the literature (Table |2) and fit these data using 
AGN and star-forming galaxy templates with the % 2 minimisation 
spectral energy distribution (SED) fitting procedure outlined in Del 
Moro et al. (2012). Briefly, we fi t the infrared data usi ng the em- 
pirical AGN template defined in Mul lanev et alj d201ll) . allowing 
for Ay =0-5 mag of extinction, with each of the "SB2"-"SB5" star- 
forming galaxy templates; the "SB1" star-formation template did 
not fit the data well in any of the sources. This provides four best- 
fitting solutions (i.e., one for each of the four star-forming galaxy 
templates) with minimum % 2 solutions. We used these solutions 
to find the star-formation rates and AGN luminosities for each of 
the sources using the methods described below. The infrared data 
points and the overall best-fit solutions are shown in Fig. [2] 

We found convincing evidence for both star formation and 
AGN activity at infrared wavelengths for all of the sources ex- 
cept for RG J0302+0010 (for which there is insufficient infrared 
data available) and SMM J0217-0503, which had no significant 
AGN component providing only an upper limit on the AGN in- 
frared luminosity. The AGN components in SMM J0943+4700 
and SMM J1237+6203 were found to be particularly bright and 
dominated the mid-infrared emission; the SED fitting results for 
SMM J0943+4700 a lso agreed well with that found from mid- 
IR spectral fitting dValiante et al.l |2007|) but no mid-IR spec- 
troscopy exists for SMM J 1237+6203. The AGN components for 
the other sources are weaker at mid-IR wavelengths. However, 
three of the sources have published mid-infrared spectroscopy 
(SMMJ0302+0010, SMM J1235 +6215, and SMM J1636+4057; 
iMenendez-Delmestre et alj|2009l) , and we found good agreement 
between the strength of the AGN component calculated from our 
SED-fitting procedure and the strength of the AGN derived from 
the mid-infrared spectroscopy. 

We used the SED-fitting results to calculate the infrared lumi- 
nosities of the star-formation components (Lvr.sf; integrated over 
8-1000 /jm) and conseq uently calculated star formation rates fol- 
lowing |Kenrucutt| 119981) ; see Table [2] This was achieved by tak- 
ing the mean luminosity of the best-fitting solutions that used the 
four star-formation templates described above. The quoted uncer- 
tainties in the star-formation luminosity and star-formation rates 
are the average of the difference between these mean values and 
the models with the highest and lowest values. We found that all 
of the sources have Lir.sf > 10 12 Lq, confirming that they are 
ULIRGs with SFRs in the range ^300-1400 M Q yr" 1 . For the 
source RGJ0302+0010 we estimated the infrared luminosity us- 
ing the radio-infrared relationship for star-forming galaxies (e.g., 



iHelou et alj[l98l) with q = 2.1 , which is a typical value for high- 
redshift ULIRGs jKovacs et al] [2006). For this source we take the 
value of Lir sp as an upper limit because the AGN will be contribut- 
ing some unknown fraction of the radio luminosity (e.g., see Fig.|2] 
which shows that some of the sources have excess radio emission 
relative to that expected from star-formation alone; also see Del 
Moro et al. 2012). 

To determine the bolometric AGN luminosities for each 
source we first calculated the AGN continuum luminosity [rest- 
frame vLv(6^m)] by taking the average AGN contribution to the 
6/jm flux of the four models described above. These average in- 
frared AGN luminosities are then converted into AGN bolomet- 
ric luminosities using the 6//m-to-2-10 keV l uminosity relati onship 
found for AGNs (conversion factor of ~0.3; lLutz et al. 20041) and 
the 2-10 keV-to-bolometric luminosity ratio in lElvis et all 
(conversion factor of ~25; see TableO. This overall correction fac- 
tor of 7.5 that we apply to convert between between vL v (6/jm) and 
AGN bolometric luminosity is in excellent agreement w ith the cor- 
rection factor of ~8±2 found bv lRichards et al.l d2006l) from their 
mean SED of SDSS quasars. We found that the two sources with 
the largest fractional uncertainties on their SED-derived AGN lu- 
minosities are also those which only have tentative evidence for 
IR-AGN activity in their mid-infrared spectra (see Table [2} • How- 
ever, we compared all the SED-derived AGN bolometric lumi- 
nosities to those predicted using the meas ured [O HI] luminosi- 
ties and the locally derived relationship of Heckman et al.l d2004h 
C^AGN=3500L[ O nj]). We found that the L[ nr] derived AGN lumi- 
nosities are on average a factor of ~ 2.5 higher than the SED de- 
rived values, with a remarkably small scatter (all sources have a 
factor of 2. 5 ±1.5), giving evidence that (1) our derived AGN bolo- 
metric luminosities are reasonable, and (2) the [O III] luminosities 
are dominated by AGN activity. To predict the AGN luminosity for 
the source RG J0302+0010, for which we were unable to perform 
the SED fitting, we therefore take the L[ 0m iderived AGN luminos- 
ity and divide by a factor of 2.5 (see Table[2}- 



4 ANALYSIS AND RESULTS 

We have studied the galaxy-integrated and spatially resolved prop- 
erties of the [O III]^A4959,5007 emission-line doublet in our sam- 
ple. Of the ten sources in our sample, six have spatially extended 
[O III] emission, two have lower quality [O III] data and limited 
spatially-resolved information and two were undetected in [O III] 
(see Table [TJ. In this section we provide the details of our analy- 
ses and give an overview of the results of the sample as a whole. 
We give specific results for individual sources in Appendix A. In 
this section we first describe the galaxy-integrated properties of the 
eight [O III] detected targets before providing a detailed analysis 
of the kinematics for the six spatially resolved sources. We defer 
discussion of the results and their implications for galaxy evolution 
to 55. 



4.1 Galaxy-integrated spectra 

In Fig. [3] we show the galaxy-integrated spectra around the 
[O III]^A.4959,5007 emission-line doublet for the eight sources 
with detected [O III] emission. The two sources that have an 
[O III] flux below our detection threshold (SMM J2217+0010 and 
SMM J22 1 7+00 1 7 ; see Appendix A9 and A 10) are not shown. The 
spectra were created by integrating over the full spatial extent of the 
observed [O III] line emission; i.e., the aperture sizes to create the 
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Figure 2. The infrared flux densities (filled circles) in the rest frame for the seven sources for which we performed SED fitting (see |J3] Tableland Del Moro 
et al. 2012 for more details). The source RG 0302+0010 (source ID number 3) does not appear due to insufficient infrared data available. For guidance, also 
shown are the overall best-fit SEDs. The total SEDs are shown as solid curves, the AGN templates are shown as dotted curves and the starburst templates are 
shown as dashed curves. The 1.4 GHz flux densities (open squares) were not included in the SED fitting process. The source RGJ0302+0010 does not appear 
due to insufficient available infrared data. These SED fits were used to constrain the contributions to the infrared luminosities from AGN activity and star 
formation activity for each of the sources (Table|2}. 



spectra were chosen by increasing the sizes until the [O III] fluxes 
reached a maximum. 

To fit the emission-line profiles we employed a % 2 minimisa- 
tion procedure, down-weighting at the positions of the strongest 
sky lines. We fit the [O III]M.4959,5007 emission-line doublet 
with Gaussian components using a fixed wavelength separation. 
The i ntensity ratio was fixed at [O IIl]X4959/[0 IIl]?i5007=0.33 
(e.g.. iDimitriievic et alj|2007l) and the width of the [O III]X4959 
line was fixed to be the same as the [O III]A.5007 line. In several 
cases the emission-line profiles are complex, with the presence of 
asymmetric wings. To characterise this complexity, we first fit a 
single Gaussian profile followed by a double Gaussian profile, only 
accepting the double Gaussian profile if it resulted in a significant 
improvement (A% 2 > 25 ~ 5a). In one source, SMM J 1636+4057, 
we find evidence for a third Gaussian component at +1350 km s 
from the narrow component; we confidently attribute this to [O III] 
as there are other emission lin es found at the same velocity as this 
component JSmail et alj|2003t also see Appendix A8). Uncertain- 
ties on the parameters are calculated by perturbing one parameter 
at a time, allowing the others to find their optimum values, until the 
X 2 value increases by 1. The mean of these upper and lower un- 
certainties are quoted throughout. The parameters for all of the fits 
and their uncertainties are given in Table [3] We note that the true 
absolute uncertainties on the flux measurements will be a factor 
of ~2-3 higher, due to unknown uncertainties on the absolute flux 
calibrations. We note that there are some discrepancies between 



the fluxes shown here and those given in Takata et al. (20o3) and 
I Alexander etal] d2010l) . We attribute this to the more reliable flux 
calibration procedures used in our work. 

In Fig. [3] it can be seen that the [O III] emission- 
line profiles are diverse across the sample, with both nar- 
row (FWHM ~ few hundred kms -1 ) and broad compo- 
nents (FWHMR^OO-lSOOkms- 1 ) which can be blueshifted 
and/or redshifted with respect to the narrow emission (—350 
< Av < 1350kms _1 ). We will discuss in Q5.l \ that the narrow 
[O III] components are likely to be tracing the host galaxy dynamics 
and merger remnants while the broad components are most likely 
due to energetic outflows (see Appendix A for a discussion on in- 
dividual objects). We will also argue that the observed diversity 
across the sample is due to a combination of orientation and obscu- 
ration effects. 

The four sources for which we were able to decompose the 
line profiles into multiple Gaussian components are those which 
have the highest flux and consequently highest signal-to-noise (Ta- 
ble |3J. It is possible that the other four detected sources also have 
both narrow and broad emission line components but we were un- 
able to decompose them in these spectra due to their lower signal- 
to-noise ratios. For example, in one source (SMMJ0217— 0503), 
we identify a broad [O III] component in a sub-region of the IFU 
datacube that is not identified in the galaxy-integrated spectrum 

In the following sections we discuss the spatially re- 
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Figure 3. Galaxy-integrated spectra, shifted to the rest-frame, around the [O III]M.4959,5007 emission-line doublet. A large variety in the [O III] emission- 
line profiles are found with both blueshifted and redshifted extremely broad (FWHMai900-1500kms _1 ) emission-line components, in addition to narrow 
emission-line components. Of particular note is SMMJ1636+4057 which exhibits two broad components separated by fsl700kms~', one blueshifted and 
one redshifted with respect to the narrow component (see Appendix A8 for details). The dashed curves show the individual Gaussian components (with an 
arbitrary flux offset) of the fits when two or three components are required. The solid curves show the best fitting overall emission-line profiles. The vertical 
dashed lines show the rest-frame wavelengths of redshifted [O III]A,4959 and [O III]^.5007 based on the redshift of the narrowest Gaussian component. The 
spectra for SMM J0302+0006 and SMM 1235+6215 are of lower signal-to-noise than the other spectra and have been binned by a factor of four for clarity. 
The parameters of the emission-line fits and their uncertainties are given in Table|3] 



solved properties of all of the [O III] detected sources except 
SMMJ0302+0006 and SMM J1235+6215. SMM J0302+0006 is 
spatially unresolved in our data while SMM J 1235+62 15 only dis- 
plays tentative evidence for extended emission (see Appendix A2 
and A6 for further details on these two sources). Due to the low 
signal-to-noise ratios of the data for these two sources it is difficult 



to determine if they are intrinsically compact or if any extended 
emission lies below our detection threshold. We do not discuss 
these two sources any further in this section. 
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Figure 4. [O III] flux, velocity and FWHM fields, from fitting single Gaussian components to individual pixels for our six sources with spatially resolved 
emission. The [O III] emission is evidently very diverse and complex with large velocity offsets (up to s=3 1700 kms -1 ) and regions of both narrow (FWHM 
» few 100 km s" 1 ) and broad (FWHM R3 700-2500 km s" 1 ) emission-line components. We attribute the narrow emission to galaxy dynamics and merger 
remnants and the broad emission to energetic outflows. Zero velocity is defined at the redshift of the narrowest component shown in the galaxy-integrated 
spectra (Fig. [3]. No corrections have been made for unknown orientation effects or dust extinction. The open circles denote the average seeing disc during the 
observations of each target and the solid bars indicate 5 kpc in extent. The dashed lines show the axes through which we produced the velocity and FWHM 
profiles shown in Fig. [6] North is up and east is left in all panels. 



4.2 Flux, velocity and FWHM maps 

In Fig. El we show flux maps, velocity fields and FWHM maps 
for the six sources for which we were able to spatially resolve 
the [O III] emission. These maps were produced using a adap- 
tive binning technique such that the high surface brightness re- 
gions are averaged over smaller regions than the outer low sur- 
face brightness regions. We do this by taking a mean spectrum 



over 3x3 spatial pixels, increasing to 5 x 5 pixels (maximum for 
SINFONI data) and ultimately to 7 x 7 pixels (maximum for NIFS 
data) if the signal was too low to produce a sufficiently high % 2 
improvement in a continuum+emission-line fit over a simple con- 
tinuum fit (A/ 2 > 16; equivalent to ~ 4o). The continuum level 
was taken to be the median of line-free continuum in the vicin- 
ity of the emission-lines. When this criterion is met, we fit the 
[O III]A,A,4959,5007 emission-line doublet with a single-Gaussian 
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Figure 4. continued. 



profile, allowing the normalisation, width and central wavelength 
to vary. We defined zero velocity at the central wavelength of the 
narrowest Gaussian component in the galaxy-integrated spectra. 

The flux maps in Fig. [4] show that the observed emission-line 
regions of these sources are diverse and irregular, with at least two 
sources showing spatially distinct emission-line regions, indicat- 
ing multiple interacting or merging systems (SMMJ0217— 0503 
and SMM J1636+4057). The [O III] emission is observed to be 
extended up to ~20kpc. The velocity fields and FWHM maps in 
Fig. [4] show that there are regions dominated by narrow [O III] 
emission (~ few hundred kms -1 ) with modest velocity gradi- 
ents (Av <200 km s~ 1 ). There are also regions which are kinemati- 
cally chaotic, extended up to 15 kpc in spatial extent, with extreme 



widths (700 < FWHM < 2500 km s" 1 ) and large velocity offsets 
(AvuptofsnOOkms- 1 ). 



The velocity fields and FWHM maps in Fig. [4] were created 
using single Gaussian components and do not separately trace the 
broad and narrow [O III] emission line components observed in 
the galaxy-integrated spectra (Fig. |3)- Instead they roughly trace 
the brightest [O III] kinematic component. For example, we will 
show that although the source SMM J 1237+6203 looks to have 
fairly quiescent kinematics in Fig. [4] the extremely broad and high- 
velocity [O III] component seen in the galaxy-integrated spectrum 
(FWHMfslOOOkms^ 1 ; Av -350km s" 1 ; Fig. is found over 
the central 5-8 kpc (see also Alexand er et al.ll2~010l) . In the follow- 
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ing sections, we therefore use further analysis methods to give a 
more complete description of the kinematics of these sources. 

4.3 Regions dominated by narrow and broad emission lines 

In Fig. [5] we show spatially-integrated spectra from regions domi- 
nated by broad or narrow emission lines. To produce these spectra 
we took our [O III] FWHM maps (Fig. [4j and integrated the spec- 
tra over the pixels that have FWHM<500kms _1 (narrow [O III]) 
and those that have FWHM>500kms _1 (broad [O III]). We do 
not do this for SMM J 1237+6203 as there are too few pixels that 
contain emission lines with FWHM>500kms _1 . For this source 
the narrow component of the [O III] emission line seen in Fig. [3] 
dominates the whole of the velocity and FWHM maps. The value 
of 500 km s -1 was motivated by the maximum line- widths gener- 
ally expected from galaxy d ynamics and merger remnants in high- 
redshift ULIRGs (see ^71 1 Alaghband-Zadeh et aT]|2012l) . We fit 
the spectra using the same method as the galaxy-integrated spectra 
(see § I4.lt and provide all of the parameters of the fits and their 
uncertainties in Table [4] We also fit any observed Hp emission in 
these regions using a single Gaussian profile and give the parame- 
ters and their uncertainties in Table |4] When no H(3 emission lines 
were detected we use the standard deviation of the emission-line 
free continuum and the FWHM of the [O III] emission-line compo- 
nents to derive 3o upper limits on the flux. 

By comparing the spatially-integrated spectra from the regions 
shown in Fig.[5]with the galaxy-integrated spectra shown in Fig. [3] 
we are able to locate the different kinematic components in the 
galaxies. For example, in the majority of the sources there is narrow 
emission (FWHMfa a few 100 kms" 1 ) extended over the full spa- 
tial extent of the emission-line regions (up to ~20 kpc). The regions 
dominated by the broadest emission lines (Fig. [5} are also the re- 
gions where the [O III] emission peaks in surface brightness and the 
[O III]/HP ratios are the largest (log([0 III]/H(3)>0.6; see Table!}, 
indica ting excitation dom inated by AGN activity in these regions 
(e.g.. lKewlevetalJl200d) . Further evidence for the broad emission 
lines being associated with AGN activity is found in the source 
SMM J02 17— 0503, which is comprised of two separate systems, 
and the broad emission line is only obser ved in the system that 
is iden tified as hosting AGN activity by Alaghband-Zad eh et al.l 
j2012l) (also see Appendix Al). In addition, SMM J 1636+4057 
displays two broad [O III] components which are co-spatial with 
the Ha broad-line region of the AGN (Menendez-Delmestre et al. 
2012; see Appendix A8). 

4.4 Velocity profiles 

To further quantify the [O III] velocity structure and spatial ex- 
tents of the different kinematic components in our sample we cre- 
ated velocity and FWHM profiles for each of the six spatially re- 
solved sources and show them in Fig. [6] To produce these pro- 
files we first repeated the pixel-by-pixel fitting routine outlined in 
ij4.2| with the addition of allowing a second Gaussian component 
to be fit at each pixel, if this resulted in a significant improvement 
(Ax 2 > 16; equivalent to Z 4o) over a single Gaussian component. 
Only two of the sources (RG J0302+0010 and SMM J1237+6203) 
had sufficiently bright narrow and broad [O III] emission-line com- 
ponents to de-compose the profiles using this method. For these 
two sources we took these double-Gaussian component FWHM 
and velocity maps and took a running median of all the pix- 
els orthogonal to a vector through the most extended emission 



(see dashed lines in Fig.|4j. For the other four spatially resolved 
sources (SMM J0217-0503; RG J0332— 2732; SMM J0943+4700 
and SMM 11636+4057) we used the single-Gaussian component 
FWHM and velocity maps shown in Fig. [4] For the sources 
SMM J02 17-0503 and RGJ0332-2732 we also plot the relative 
velocities of the Ha emission, extracted from the regions shown in 
Fig. [5] on these velocity profiles. For the source SMM J1636+4057 
we plot the [O III] FWHM and velocity offset values from the 
spectra extracted from the broad [O III] region (Fig. |5} and 
Ha and CO(l-O) data from the literature dSwinbank et aLH2005t 
llvison et"al]|201 ll ; Menendez-Delmestre et al. 2012; see Appendix 
A8 for more details). 

Fig.[6]enables us to separately trace the kinematics of the nar- 
row (FWHM a a few hundred km s" 1 ) and broad (FWHM si700- 
1400 kms~') [O III] emission-line components. The narrow emis- 
sion lines have small velocity gradients (Av <200kms _1 ) and ap- 
pear to be at roughly the same velocity as the Ha emission and 
molecular gas, when these data are available for comparison. In 
contrast, the broad components are offset from the narrow emis- 
sion by A|v| ~ 150-1400 km s and are found over ~4-15kpc 
in observed linear extent (except SMM J0217— 0503 for which 
the broad emission is spatially unresolved; ^4 kpc). The source 
RGJ0332— 2732 has emission that is not clearly separated into 
broad and narrow components; however, we note that the emis- 
sion is broadest (FWHM~700kms~') in the brightest central re- 
gion. The source SMM J0217— 0503 consists of two separate sys- 
tems, which are indicated in Fig. [6] In the northern system we do 
not observe any narrow emission lines; however, the [O III] emis- 
sion has a large velocity offset from the peak in the Ha emission 
(Av = —460 ± 60 km s ; see Appendix Al for further details). Of 
particular note is the source SMM 1636+4057 for which the broad 
components are found at v ~ ±850kms~' from t he systemic ve- 
locity of the the Ha broad-line region of the AGN dSwinbank et al.l 
2005; Menendez-Delmestre et al. 2012; see Appendix A8). 

Due to the large uncertainties in the velocity offsets and 
FWHM of the broad [O III] components, we do not interpret their 
kinematic structure any further than quoting their spatial extent and 
velocity offsets from the systemic redshifts. However, we did per- 
form tests to ensure that the broad emission lines are truly spatially 
extended and that the observed extents are not just the result of in- 
trinsically compact regi ons being extended by th e seeing during our 
observations. Following I Alexander etafl | |2010|) we did this by ex- 
tracting spectra from different regions and we found velocity struc- 
ture in the broad emission lines, confirming that they are spatially 
extended (see Appendix A for details on individual sources). 

For the remainder of this work we quote the spatial extent of 
the broad [O III] emission as the extent of the broad [O III] regions 
shown in Fig. [5] (see Table f4j. These values are consistent, within 
the quoted uncertainties, with spatial extents of the broad [O III] 
emission found in the velocity profiles shown in Fig|6] The source 
SMM J 1237+6203 does not appear in Fig.[5](see jg3) , we therefore 
take the spatial extent from Fig. HQ (~5 kpc) which is conserva tive 
based on the 4-8 kpc extent found by lAlexander et ajj ( 120101) for 
this source. The quoted spatial extents are likely to be lower limits 
as they are not corrected for orientation and low surface bright- 
ness broad [O III] emission will not be detected by our fitting pro- 
cedures (described in ^4/2}. We define the velocity offsets of the 
broad emission lines from the systemic, as the velocity differences 
between the broad and narrow [O III] emission lines in the spectra 
extracted from the broad [O III] regions (see Fig. [5] and Table [4}. 
For SMMJ0217-0503 and RGJ0332-2732 we use the velocity 
offset to the narrow Ha emission, as no narrow [O III] emission 
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Figure 5. Left: Maps of our spatially resolved sources colour-coded to separate regions dominated by narrow [O III] emission lines (green: FWHM 
< 500kms~') or broad [O III] emission-lines (orange: FWHM > 500 kms~') determined from the FWHM maps shown in Fig. [4] For the source 
SMM J 1237+6203 the narrow component dominates the [O III] flux throughout and therefore we are unable to define the separate spatial regions using 
this method (see i]4.3l for details). The narrow emission lines are observed over the largest extents but the broad emission lines are still observed in extended 
regions («4-15 kpc). The open circles denote the average seeing discs during the observations, the solid bars indicate 5 kpc in extent and the crosses indicate 
the peak in [O III] flux. North is up and east is left. Middle and Right: The spatially-integrated spectra around the [O III]M.4959,5007 emission-line doublet 
extracted from the regions indicated in the left panels. The central regions exhibit complex emission-line profiles. The dashed curves show the individual 
Gaussian components (with an arbitrary flux shift) of the fits when two components were required. The solid curves show the best fitting overall emission-line 
profiles; the parameters and their uncertainties are shown in Table|4] The vertical dashed lines shows the centroids of the [O III]M4959,5007 narrow line 
components in the galaxy-integrated spectra (Fig. [3)- The vertical dotted lines for SMM J0217— 0503 and RGJ0332— 2732 show the expected position of 
[O III] given the redshift of the Ha emission-line from the same spatial region (Alaghband-Zadeh priv. comm.). Fo r the source SMM J163 6+4057 the vertical 
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Figure 6. A velocity and FWHM profile of the [O III] emission lines for our six spatially resolved sources, along vectors at the angles indicated in the 
upper-right of each panel (see £|4.4| for how these were produced). When it was possible to de-couple the [O III] emission-line profiles into two components 
we have shown these as filled circles (narrow) and open squares (broa d). In some cases Hcc data and CO(l-O) data are plotted, taken from the literature 
fSwinbank etail2005t Ivison et al. 2011; Alaghband-Zadeh et al. 2012; Menendez-Delmestre et al. 2012). The data points for the broad [O III] components 
for the source SMMJ1636+4057 are taken from the fits to the spectra shown in Fig. [5] In all cases there are underlying narrow [O III] emission lines (FWHM fs 
few hundred kms~') with small velocity gradients (Av <200kms~') which we associate with galaxy dynamics and merger remnants. Broad [O III] emission 
lines (FWHMKS700-1400 km s~ 1 ) are offset from these galaxy kinematics which we interpret as energetic outflows. The error bars indicate are representative 
lcj uncertainties for the [O III] measurements. The horizontal solid lines indicate the approximate uncertainties in the spatial positions of the measurements 
when applicable. The dotted lines indicate zero velocity using the redshift of the narrowest [O III] components in the galaxy-integrated spectra (Fig. [5). 
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Figure 7. The [O III] FWHM from the galaxy-integrated spectra of our eight 
[O III] detected targets and high-redshift ULIRGs from the literature versus 
1.4 GHz radio luminosity. Both broad components (FWHM>700kms _1 ; 
dashed line) and narrow components of the profile fits have been plotted (see 
Table[3}. For SMMJ0217— 0503 we have plotted the northern and south- 
ern systems separately (see Table [4}- We also plot the galaxy integrated 
Ha FWHMs from other high-redshift ULIRGs with simil ar radio luminosi- 
ties JSwinbanket~aL. 2004; Alaghband-Zade h et all20l3) . The sources that 
are identified as AGN (using methods not based on emission-lines widths), 
have been highlighted (see ij lS.ll for details). The non-AGN ULIRGs have 
FWHMs that lie below w700bns _1 . The narrow [O III] components of 
our sample are consistent with the emission-lines widths from non-AGN 
ULIRGs and presumably trace galaxy dynamics and mergers while the 
broadest components appear to be constrained to the AGN host galaxies 
and are likely to be tracing AGN-driven outflows. 



is observed. For SMM 1636+4057 we use half of the velocity off- 
set between the two broad components, which is also consistent 
with the offset between each component and the velocity of the Ha 
broad-line region (see Appendix A8). 



5 DISCUSSION 

Our IFU observations of eight z = 1.4 - 3.4 ULIRGs that host AGN 
activity have revealed extreme gas kinematics with very broad 
[O III] emission (FWHMf»700-1400kms _1 ), extended up to 
15 kpc and with high-velocity offsets, up to ~850 km s~ 1 . We have 
also identified narrow [O III] emission lines with moderate widths 
and velocity gradients (FWHM<500kms _1 ; Av <200kms _1 ) 
over the full extent of the observed emission-line regions (up to 
~ 20 kpc). In the following discussion, we focus on the six sources 
for which we have high-quality spatially-resolved data and show 
that the narrow emission lines are consistent with tracing galaxy- 
dynamics and mergers, while the broadest emission lines are con- 
sistent with being predominantly due to energetic outflows. We 
then go on to discuss the impact that these outflows could be having 
on the future growth of the central BHs and host galaxies and the 
implications of these results for the formation of massive galaxies. 

5.1 Tracing galaxy dynamics, mergers and outflows 

In six of our sources we find an extended [O III] emission- 
line region (up to ~20kpc; Fig. |4), which are often found in 



gas-rich systems containing luminous AGN (e.g., IColina et al.l 
ll999tlNesvadba et alj2o61lGreene et al.l201ll : Villar-Martm et al. 
201 la; lMatsuokj|2012h . The presence of these emission-line re- 
gions enables us to trace the kinematics of the ionised gas in these 
galaxies. 

We identified narrow [O III] emission (FWHM up to a few 
xlOOkms -1 ) with irregular velocity profiles and small velocity 
gradients (Av < 200kms~ I ) over the full extent of the emission- 
line regions (Fig.[4]and Fig.[6). The line widths and velocity offsets 
observed in this narrow [O III] emission are simil ar to the kinemat - 
ics of gas found in simulated merger remnants dCox et al.| [2006). 
The kinematics of these narrow emission lines are also broadly con- 
sistent with other high-redshift star-forming galaxies that have sim- 
ilar radio-lu minosities, man y of which are also underg oing mergers 
(e.g. ISwinbanket al.ll2004 Nesvadb a et al. 2007b; iLehnert et al.l 
120091; I Alaghband-Zadeh et al.l l2012h as well as nearby ULIRGs 
(e.g., IColina et alj|2005l) . To illustrate this, in Fig. [7j we compare 
the galaxy-integrated emission-line FWHMs of our AGN sample 
to the Ha emission l i ne FW HM in other high-redshift ULIRG s 
from ISwinbank et"aH d2004l) and lAlaghband-Zadeh et al.l d2012h . 
We highlight sources that are identified as AGN on the basis of their 
emission-line ratios, X-ray emission, UV spectral features or their 
mid-infrared spectra. Fig. [JJ shows that the narrow [O III] emis- 
sion lines (FWHM^SOOkms- 1 ) in our sources are comparable to 
the emission-line widths of non-AGN high-redshift ULIRGs and 
therefore are likely to be tracing the same kinematic components 
seen in these sources (i.e., galaxy dynamics and mergers). The ve- 
locity structure of t he emission lines seen in the IFU data of the 
non-AGN ULIRGs 1 Alaghband-Z adeh et "ai]|2012l) are also similar 
to the narrow [O III] components observed here. In particular, for 
the two sources where there are IFU observations covering both the 
[O III] and Ha emi ssion lines (SMM J02 17-0503 and RGJ0332- 
2732; this work and lAlaghband-Zadeh et al.|[2012h the kinematics 
of the narrow [O III] and Ha emission-line components are roughly 
consistent (e.g., see Fig. [5}. 

In contrast to the narrow emission lines, the identification 
of broad emission lines (FWHM^700kms~') in the galaxy inte- 
grated spectra of high-redshift ULIRGs appears to be constrained 
to galaxies that contain known AGN activity (Fig. [7j|3 This is in 
agreement with studies of low-redshift ULIRGs where the broadest 
emissi on lines are found in systems hosting powerful AGN activity 
(e.g. IVeilleux et aljfl999l IZheng et alj|2002l ; IWestmoquette et all 
120121) . Although broad emission-line components may be found 
in non-AGN high-redshift star-forming galaxies, they are gener- 
ally faint and are difficult to identify other than by stacking the 
spectra of multiple s ources or with excelle nt quality data (e.g ., 
ISwinbank et al.l2004l ; Nesvadba et al. 2007b; lLe Tiran et al.l201lh . 
However, we note that we are basing these arguments on small 
samples of galaxies which are biased in their selection and to 
fully address if bright and broad emission lines, such as those in 
our sources, are unique to galaxies hosting AGN activity, requires 
larger systematic samples of spectra and IFU data. In addition to 
comparisons with the literature, there is further evidence that the 
broadest emission lines are associated with AGN activity in our 
sample. The broadest emission lines are found to be spatially co- 



3 The one source from ISwinbank et all |2004|) with FWHM>700kms"' 
that is not identified as an AGN (Fig. |7j has a complex spectrum and has 
not been covered by multi-wavelength observations in the literature. It is 
therefore not possible to rule out the presence of a significant AGN in this 
galaxy. 
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incident with the most likely location of the AGN in these galax- 
ies based on very high [O III]/HP emission-line ratios (see E|4.3b . 
The most striking example of this is the source SMM J1636+4057 
where the broadest [O III] emission lines are found to be co-spatial 
with the AGN, and not in the region of unobscured star formation 
(details in Appendix A8). 

Broad and asymmetric [O III]XX4959,5007 emission-line pro- 
files, such as those found in this work, are most commonly 
attributed to high-velocity outflowing gas (e 



198ll:lGreene & Holl2005l:[Barth et al.f 2008; Nes vadba et a"Hl2008l ; 



Greene et alj 201 lh . In particular IFU observations of local 



ULIRGs that host AGN activity reveal broad and spatially extended 
emission lines (simil ar to those that we observe) associated with 
energ etic outflows l lRupke & Veilleuxl 1201 lh IWestmoquette et all 
l2012h . On the basis of the above arguments the broadest kine- 
matic components (FWHM>700kms _I ) appear to be associated 
with AGN activity and therefore do not appear to be the result 
of mergers and galaxy dynamics or even (as we will further ar- 
gue below) star-formation driven outflows, and therefore AGN- 
driven outflows is the most natural explanation. We find con- 
vincing evidence for extended (> 4kpc) outflows in four of our 
sample (RG J0302+0010; SMM J0943+4700; SMM J 1237+6203 
and SMM J 1636+4057). The source RG J0332-2732 only has very 
broad [O III] emission lines in the most central regions, and the 
broad [O III] emission line in the source SMM J02 17-0503 has 
a low signal-to-noise ratio and is spatially unresolved. These two 
sources have lower quality data than the other four spatially ex- 
tended sources and the evidence for outflows in these objects is 
weaker (see Appendix Al and A4 for further discussion on these 
two sources). 



5.2 Outflow properties 

We can use the measured properties of the broad [O III] emission to 
place constraints on the outflow properties. In Fig.|8]we plot the ob- 
served spatial extent of the broad [O III] emission lines against the 
velocity offset from the systemic redshifts of the host galaxies (val- 
ues are from the fits to the spectra shown in Fig.[5]and are provided 
in Table[4]l0 The dashed lines in Fig.[8]represent estimated energy 
injection rates which we discuss later. The broad [O III] emission- 
line properties (Av ~50-850km s ; spatial extent ~4— 15 kpc) can 
be explained by the presence of an AGN-driven bi-polar outflow 
plus the addition of obscuring material in the host galaxy. In Fig. [9] 
we show a simple schematic diagram to attempt to explain the 
properties and diversity of the broad [O III] emission lines. Hydro- 
dynamical simulations and analytical models predict that a pow- 
erful AGN-wind, launched from the central BH, cou ld sweep up 
clouds of gas resulting in an kpc-scale outflow (e.g.. lDebuhr et al] 
l2012l ; IZubovas & Kind[20ll and references therein). Galaxy-wide 
outflows are thought to be expanding bubbles on either side of 
the host galaxy, being forced into this sha pe by the denser mate- 
rial in the centre of the host-galax y (e.g., IZubovas & Kindl2012l ; 
Faucher-Giguere & Ouataerdl2012l) . As we are observing the out- 
flow through luminous [O III] emission, the observed properties of 
the outflow will also depen d on the size and orientation of the AGN 
ionisation cones (e.g. JCrenshaw e t al. 2010). Multiple clouds being 



4 The source SMM J 1237+6203 does not appear in this table. For this 
source we use a spatial extent of 5 kpc (see Fig. [6j and the FWHM and 
velocity offset from the galaxy-integrated spectrum (Table [3). See the end 
of q4.4l for details. 



caught up in the bulk flow of the outflow (e.g., IZubovas & Kind 
120121) or perturbed ionised gas would result in broad [O III] emis- 
sion lines (e.g., Villar-Martm et al. 2011b). The 'near-side' of the 
outflow would therefore be observed as a broad and blueshifted 
emission line and the 'far-side' would be observed as broad and 
redshifted emission (see Fig. [9}. However, part of the outflow could 
be obscured by material in the host galaxy, as is likely to be the case 
in clumpy and du s ty systems such as h i gh-redshift ULIRG s (e.g., 
Ilvison et al.ll200l iTacconi et al.l [200^ : Engel et ai]|20ld) . High- 
redshift ULIRGs are kinematically disturbed, potentially resulting 
in a complex distribution of obscuring material, which, along with 
the relative orientation of the AGN, could explain why in some 
cases we only see either a blueshifted (i.e., when only the near-side 
of the outflow is observed) or redshifted broad emission line (i.e., 
when o nly the far- side of th e outflow is observed; also see Fig. 6 
and 7 in ICrenshaw et aill2010l) . 

The range in spatial extents and velocity offsets of the implied 
outflows that we observe (Fig. [8} can also be explained using this 
simple model (Fig.[9]l. The observed velocity offsets and observed 
spatial extents will be highly dependent on the orientation of the 
outflow with respect to the line of sight. If the axis of the outflow is 
orientated along the line of sight, high velocity offsets and a small 
spatial extent would be observed (e.g., this could explain the prop- 
erties of SMM J1636+4057). Conversely, if the axis of the outflow 
is in the plane of the sky, a small velocity offset and a large spatial 
extent would be observed (e.g., this could explain the properties of 
SMM J0943+4700); see Appendix A for further discussion on indi- 
vidual sources. In contrast, the observed FWHM will be less depen- 
dent on the orientation. We emphasise that this is a very simplified 
model and other interpretations could also be valid. For example, 
the observed emission-line properties will also be dependent on the 
intrinsic properties of the outflow, such as their intrinsic velocity 
and distance traveled from the BH. 

We have argued that the measured velocity offsets observed 
may not represent the true velocity of the outflow due to an ori- 
entation effect. We attempt to estimate corrected velocity offsets 
by taking the two sources which contain a known broad-line AGN 
(SMM J 1237+6203 and SMM J1636+4057; see Appendix A for 
details). In these sources the identification of the broad-line AGN 
component suggests that the outflow could be primarily orientated 
towards us and therefore the observed velocity offset is likely to 
be close to the true outflow velocity. For these two sources we find 
that the average of FWHM/Av is ~ 2; reassuringly we also find a 
median ratio of FWHM/Ai' ~ 2 for low-redshift broad-line AGN 
that display [O III] properties indicative of outflows (data from 
Mullaney et al. in prep). We therefore suggest that FWHM/2 may 
be an adequate approximation of the true outflow rate, although 
subject to large uncertainties. Using this as a diagnostic, we es- 
timate bulk outflow velocities in the range of ~300-700kms _1 , 
which are similar to the predicted kpc-scale outflow velocities pre- 
dicted by hydrodynamical simulations and analytical models of 
outflows (~300-1500kms _1 ), driven by AGN in objects with sim- 
ilar properties to the sources in our sample (i.e., ULIRGs, with 
/-aon ~ 10 46 erg s~ 1 ; iDebuhr et al J |2012| ; IZubovas & Kind |2012| ; 
iFaucher-Giguere & Ouataerul2012l) . 



5.2.1 Estimating outflow energy rates and determining what is 
powering the outflows 

To attempt to quantify the impact the outflows may be having on 
the host galaxies we consider approaches to estimate the energy 
injection rates into the ISM. We do not have good constraints on 
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Figure 8. The projected spatial-extent of the broad [O III] emission-line regions against velocity offset from the systemic redshift (left) and FWHM (right) for 
our spatially resolved sources (see Table |5]and Table |4j- The numbers inside the symbols correspond to the source IDs. Also shown are other high-redshift 
sources with IFU data (symbols are the same as in Fig.[T]. If two sides of an outflow are observed, half of the velocity difference between them is used, and the 
average FWHM of each side. For spatially unresolved data we use the seeing disk as an upper limit of the extent. For illustrative purposes, the dashed lines show 
the expected energy injection rates (E) using a simple outflow model (see £|5.2| for details). Strong evidence for extended outflows is not confirmed for source 
IDs 1 and 4 (see t|5.U . Despite 3-4 magnitudes difference in radio luminosities, HzRGs and high-redshift ULIRGs containing AGN activity have comparable 
velocity offsets and emission-line widths and are both potentially dumping energy into their host galaxies at considerable rates (calO 43 — 10 46 erg s ). 



mass of the gas entrained in the outflows or on the geometry of the 
outflow, making such estimates uncertain. For illustrative puiposes, 
we consider the energy injection rates assu ming an energy con- 
serving bubble into a uniform medium (e.g . , iHeckman et al.lll990l ; 
IVeilleuxetaLll2005l :l Nesvadba et al.|[2008i : and references therein) 
which gives the following relation: 

1.5 x 10 4 %v3 ()00 n .5 ergs- 1 , (1) 

where no is taken to be half the extent of the observed broad [O III] 
emission (in units of lOkpc) and Viooo are the velocity offsets in 
units of 1000 km s -1 (see Table[4]l. The ambient density (ahead of 
the expanding bubble), M0.5, is in units of 0.5 cm -3 for which we 
have assumed rco.5 = 10 

In the left panel of Fig. [8] we plot on tracks indicating fixed 
energy injection rates using Equation Q] We have already dis- 



5 This density assump ti on is based on the the valu e s of n p 5=l^t used 
in e.g. IHeckman etal] jl990) and iNesvadba et~all (2008). This value 
can be indirectly est imated from the (post-shock) electron density (see 
IHeckman etal] 1 19901 and references therein for further details), derived 
from the [S IIjM.6171,6731 emission-line doublet which is sensitive to 
electron densities in the range 100-10 5 cm -3 IQsterbroc k 1989). We used a 
stack of the res t -frame optical spectra of the high-redshift ULIRGs from 
ISwinbank et alj J2004I) to derive a flux ratio o f Rpn] = fi/)7ifi/ffii/iT ii = 
1.1 ± 0.3. This is consistent with the HzRGs of Nesvadba et al. 200j| who 
find R[sn]=10-l.l, such that the implied densities of HzRGs and high- 
redshi ft ULIRGs are simi lar. We therefore use the same value of no. 5 = 1 
as in INesvadba et al. 12008). Although we have no direct way of measuring 
no. 5 this unknown will only vary our order-of-magnitude estimates of the 
energy injection rates by a factor of a few. 



cussed that the observed velocity offsets are unlikely to be rep- 
resentative of the intrinsic velocities and that FWHM/2 is plau- 
sibly a better estimate of the true outflow velocity. Therefore, 
in the right panel of Fig. [8] we plot the extent of the broad 
[O III] emission against FWHM and tracks of constant energy 
injection rates using Equation [T] but replacing the velocity off- 
set with FWHM/2. For the two sources with a known broad-line 
AGN in this sample (SMM J 1237+6203 and SMM J1636+4057) 
and the sources with broad-line AGN in the co mparison samples 
(i.e. lNesvadba et al hood ie ano-Diaz et al .120121) . the two diagnos- 
tics for the bulk outflow velocity give energy injection rates that 
are consistent within a factor of a few. We establish that the spa- 
tially extended outflows in these high-redshift ULIRGs are poten- 
tially injecting energy into their host galaxies at considerable rates 
of ss(0.08-3)x 10 45 ergs" 1 , where the range indicates the values 
for different sources . Over a typical AGN lifetime of 30Myr (e.g., 
Hopki ns et al. I l2005h the total energy injection would be of order 
~(0.8-28)x 10 59 ergs; however, we note that outflows could con- 
tinue ~10 times longer than the active time of the BH dKing et aT] 
l201lh which would increase this energy injection by another order 
of magnitude. These outflow energy injection rates are comparable 
to the typical binding energy of the host galaxies in these systems 
of aslO 59 ergs (assuming a spheroid mass of Ri lO 11 M^, a velocity 
dispe rsion o ~200kms-' and i? e =4kpc; e.g. jBinnev & Tremaind 
Il987l) and demonstrate that these outflows could unbind at least a 
fraction of the ISM from the host galaxies. 

We next consider what could be powering the outflows based 
on energetic arguments. In Fig.[l0]we compare the ratio of the out- 
flow energy rates with the potential energy input rates from: (1) 
the radiative output of the AGN (i.e. the AGN bolometric luminos- 
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Figure 9. A schematic diagram to illustrate a possible interpretation of 
the observations of the broad [O III] emission lines. If a fast wind is ini- 
tially launched around the central BH it could sweep up clouds of gas in 
a bi-polar outflow (as assumed by some models e.g., Debu hr et al.ll2012l ; 
IZubovas & Kinj20ll see § l5.2l for more discussion). The orientation with 
respect to the line of sight will determine the velocity offsets and spatial 
extent observed; however, broad emission lines will be always be observed. 
This could explain the diversity in the properties of our observations. In 
addition, obscuring material in the host galaxy would result in parts of the 
outflow not being visible which could explain why in some sources we only 
observe either blueshifted or redshifted broad emission lines. 



ity; Z-agn) an d (2) energy output from star formation (£sf)- We 
calculate outflow energy injection rates using Equation 1, replac- 
ing v with FWHM/2 and r as half the spatial extent of the broad 
[O III] region (see Table [4] also see footnote |4). We have esti- 
mated the bolometric luminosity of the AGN using SED modeling 
(described in ij3j values given in Table [2}. We estimate the poten- 
tial energy input from star-formation using the star formation rates 
from our SED modelling (Table[2) and assuming that the mechani- 
cal energy injection from supernovae and stellar wi nds would be at 
most mi x 10 41 (SFR/MfT) yr~ 1 ) erg s~ 1 (following iLeitherer etail 
ll999l and lVeilleux et al]2005l) FlWe note that the assumptions made 
are highly idealistic, and we have not, for example , considered the 
case o f a mo mentum conserving outflow. However. iNesvadba et al.l 
< l2006t 2008) follow several methods to estimate the energy in- 
jection in similar outflows in HzRGs and find that the different 
methods gen erally agree on the order of magnitude level (also see 
IVeilleux et al1l2005l for further discussion on the energetics of out- 
flows and the potential energy input from star formation and AGN 
activity). 

Fig. 1 101 demonstrates that the radiative output of the AGN is 
energetically viable as the dominant power source of the outflows 



6 If we were to follow instead Dalla Vecchia & Schavd 120081) to estimate 
this energy input, the values would be a factor of «2 lower. Conversely, 
radiation pressure from star-formation could potentially contribute a com- 
parable amount of p ressure to stellar and supernovae winds in ULIRGs 
jveilleux et alj|200a) . As we are considering order of magnitude estimates 
only, these uncertainties do not affect our conclusions. 



observed in all of these sources. Assuming the bolometric luminos- 
ity provides an estimate of the initial energy input from the AGN, 
so long as ~0.2-5% of this energy can couple to the ISM, the radia- 
tive power of the AGN in these systems is large enough to drive the 
observed outflows. This range of coupling efficiencies is similar to 
the values required by many models to reproduce the properties of 
local massive galaxies (e .g., ISpringel'e t al. 2005; Di Matteo et al.l 
l2005l : iDebuhr et al.ll2012h . In contrast, the required coupling effi- 
ciencies for star formation are an order of magnitude higher (~ 1 5% 
to 2>100%). For the three sources with implied coupling efficien- 
cies of ~ 15-50% star formation may play a significant role; how- 
ever, based on the observed properties of the broad [O III] emission 
presented earlier (see i|5.1b . and the energetic arguments presented 
here, we tentatively suggest that the AGN may play a dominant role 
in powering the outflows. Given the current quality of the available 
data, a more detailed analysis is not yet warranted and the values of 
input energetics should be used with care. 



5.2.2 The potential fate of the outflowing gas 

Based on energetic arguments we have already demonstrated that 
the outflows we observe are likely to unbind the gas from the 
host galaxies. However, many models of galaxy formation re- 
quire powerful AGN-driven outflows to remove the gas completely 
from the host galaxy potentials to efficiently truncate star forma- 
tion and remove the fuel for future accre tion onto the BHs (e.g., 
|Pi Matteo et"alj 2005; Ho pkins et aiT2 008). Will the outflowing gas 
we observe be able to escape the host galaxy and consequently re- 
move the gas supply for further star formation and future accretion? 
In Fig. [10] we plot the broad [O III] velocity offset against FWHM 
for our spatially resolved sources and also plot other high-redshift 
AGN with available IFU data from the literature. By taking an es- 
timated galaxy mass of lO^Mfp (typical of these source s; e.g., 



1p (typi* 

lEngel et alj2"oiol : IWardlow et al.l201 ltlHainline et al.l201 ll) the es- 
cape velocity at a distance of 4kpc would be v e .s C ~460kms 
(see Fig. |10t . Considering the observed velocity offsets (Av=50- 
850kms _1 ), only the outflows in two sources (SMM J 1636+4057 
and SMM J0217— 0503) appear likely to exceed this escape veloc- 
ity. However, if we instead consider the the maximum velocities of 
the emission-line g as (Vmax 

-400-1400 km s" 1 ; defined 

as v max — 

|v|+FWHM/2; i.e., iRupke et alj|2005al fbh then the outflow veloc- 
ities are found to exceed the escape velocity (see dashed curve 
in Fig. I lOt in all sources except RG J0332— 2732, suggesting that 
some fraction of the gas is able to escape the host galaxy potential. 

Unless the gas is completely removed from the host galaxy 
potential, or is unable to cool on short timescales, significant star 
formation and BH accretion may be able to recommence after the 
outflow episode. We therefore consider a simplified approach to in- 
vestigate if the outflowing gas will be able to escape the dark-matter 
halo. Assuming a static dark matter halo with a m ass of ^10 13 M Q 
(expected for thes e sources; Hickox et al. I l2012l) . a Navarro Frenk 
and White (NFW; lNavarro et"aL"H "l996) density profile, and a point 
mass of ?sl0 Mq (i.e., the host galaxy), a particle initially at 
~4 kpc from the centre of the host galaxy would need to be trav- 
elling at v ~ 1000 km s -1 to escape beyond the virial radius of the 
halo (~440kpc; M. Lovell priv. comm.; see Fig. UOt. On the ba- 
sis of this simple model, where we assume the effects of gravity 
only, particles travelling at ~400-800kms~' at ~4kpc will return 
to the host galaxy within a few hundred Myrs. From this analy- 
sis only two sources have maximum outflow velocities that could 
potentially unbind some fraction of the gas from the dark mat- 
ter halo (Fig. Hot , which suggests that the gas in the majority of 
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Figure 10. Left: The ratio of the estimated outflow energy rate (£ ut) to the derived input energy rate from AGN (£agn) against the ratio of the outflow energy 
rate to the predicted energy input rate from star formation (£"sf) for five of the six sources with spatially resolved [O III] emission. For the sixth source we 
only have upper limits for both £agn and £ i,t and therefore have no constraints on the ratio of these quantities. The numbers inside the symbols correspond 
to the source IDs. We note that we lack strong evidence for extended outflows for source IDs 1 and 4 (see gj]) The shaded regions indicate where a >100% 
coupling efficiency is required between the input energy and the gas to power the outflows, effectively ruling out such an input energy mechanism. The 
bolometric luminosity of the AGN appears to be energetically favourable to power all the outflows requiring only «0.2-5% of the bolometric output to couple 
to the gas; however star formation cannot be ruled out in all cases. It should be noted that the values used rely heavily on assumptions and should be taken 
as illustrative only, with uncertainties at the order of magnitude level. Details of the derived values that are plotted are provided in i]5.2l (also see Table[2j. 
Right: Broad [O III] emission-line velocity offset versus FWHM for the spatially resolved sources in our sample (values taken from regions of broad [O III]; 
see Table|4j. Also shown are other high-redshift sources with IFU data as described in Fig. [T] The maximum velocities (v raax = |v| + 4 FWHM) of the gas 
are potentially high enough to exceed the typical escape velocity of a typical massive galaxy (10 Mq at 4kpc; dashed curve). However, using simplistic 
arguments (see £|5 .2t only one of our sources has a maximum outflow velocity that could potentially unbind some of the gas from a typical galaxy halo 
(»10 13 Mq at «4kpc; dotted curve). 



the systems may eventually cool and re-ignite another episode of 
AGN activity and star formation. However, it has been proposed 
that quasar mode outflows, such as those observed here, may be 
an effective way of "pre-heating" the gas at early times prior to 
a mechanical-heating dominated phase (i.e., the so-called radio 
mod e), which will then efficiently prevent the gas from cooling 
(e.g jvan de Voort et al.l201 lMMcCarthv et alj201 ll) . The outflows 
we observe could be a means of energetically pre-heating and eject- 
ing gas at high redshift and therefore could be a crucial stage in the 
evolution of massive galaxies. In this scenario these galaxies will 
therefore be the progenitors of the radiatively weak, low-redshift 
AGN in massive galaxies which are obser ved to have mechanical 
heating that is preventing gas from cooling dBest et al.l2005l [20061 ; 
ISmolcic et ai]|2009l : iDanielson et alj2 012). 



6 CONCLUSIONS 

We have taken Gemini North-NTPS and VLT SINFONI IFU ob- 
servations around the emission-line doublet [O III]A4959, 5007 of 
eight high-redshift ULIRGs that host AGN activity. Our aim was to 
search for broad, high-velocity and spatially extended [O III] emis- 
sion as a tracer of powerful AGN-driven outflows over galaxy- wide 
scales that are predicted by the quasar-mode in galaxy evolution 
models. We have strong evidence that such outflows exist in at least 
a fraction of the high-redshift ULIRGs. Our main conclusions are 
the following: 

(i) We find narrow [O III] emission lines (FWHM~ few hun- 
dred kms -1 ) extended over galaxy- wide scales (up to ~20kpc). 



By comparing these observations with high-redshift ULIRGs that 
are also undergoing intense starbursts and mergers but do not host 
significant AGN activity we show that the narrow [O III] emission 
lines we observe (FWHM up to a few lOOkms -1 ) are consistent 
with tracing host galaxy dynamics and merger remnants (i j5.lt . 

(ii) In the four most luminous sources we find extremely broad 
[O III]M4959,5007 emission (FWHMa^OO-WOOkms- 1 ) with 
high velocity-offsets (up to ~850kms _1 ), extended over «4— 
15kpc; i.e., the key signatures of galaxy-scale energetic outflows 
( ET) . The other four sources exhibit broad [O III] emission but 
we have limited constraints on the spatial extent of the broad [O III] 
emission due to lower quality data. 

(iii) We apply a simple outflow model to show that the spatially 
extended outflows are potentially moving through the ISM with an 
energy rate of ~(0.08-3)x 10 45 erg s _1 , sufficient to unbind some 
of the gas from the host galaxy and potentially disrupting future 
accretion and star-formation in the host galaxy ( g|5.2.1| (. These or- 
der of magnitude estimates are model dependent and rely on esti- 
mates of the density and therefore should be considered as illustra- 
tive only. 

(iv) By estimating the energy injection rates from the AGN and 
star formation driven winds we find that the outflows require ~0.2- 
5% of the AGN luminosity to couple to the gas, while star forma- 
tion require coupling efficiencies of 15%-3>100% ( ^5.2.11 1. We use 
a combination of energetic arguments and comparison to non- AGN 
systems to show that, although star formation may also play a sig- 
nificant role in some sources, the AGN activity is feasibly the dom- 
inant power source for driving the outflows in all sources. 

(v) The high maximum outflow velocities observed in these 
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galaxies suggest that at least some fraction of the gas will be un- 
bound from the host galaxies, but significant fractions of gas are 
very unlikely to be completely unbound from the galaxy haloes 

(EH). 

The galaxies observed here are two to three orders of magnitude 
more common than the extreme high-redshift radio galaxies that 
are often associated with powerful AGN large-scale outflows at 
high redshift. We have shown that energetic large-scale outflows 
are present in at least a fraction of high-redshift ULIRGs that host 
AGN activity, which are believed to be the progenitors of today's 
massive elliptical galaxies. We have suggested that these outflows 
may be a crucial phase in the evolution of massive galaxies. 
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Target SED Data and Derived Properties 



ID 


Source 


s t 


S24 


S70 




S350 


S500 


^850 




SFR 


^AGN.SED 






C"Jy) 


(mjy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(1O' 2 L ) 


(M yr-') 


(10 46 ergs-' 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(ID 


(12) 


1 


SMM J0217-0503 


44.7±2.3 


0.49±0.05 




36.3±5.8 


33.9±6.3 


21.7±6.8 


4.4±1.7 


4.5±0.8 


780±160 


<0.2 


2 


SMM J0302+0006 


63.4±6.8 


0.45±0.05 


<13.6 








4.4±1.3 


2.0±0.3 


340±50 


0.14±0.09 


3 


RGJ0302+0010 




<0.24 










<4.5 


<2.1J 


<360 


2.0±0.2** 


4 


RGJ0332-2732 


7.8±0.2 


0.94±0.09 




27.5±5.8 


27.5±6.3 


20.3±6.8 


<3.6 


4.9±0.9 


850±120 


2.2±0.4 


5 


SMM J0943+4700* 




0.57±0.06 




22.6±5.8 


26.5±6.3 


17.3±6.8 


8.7±1.5 


7.5±1.5 


1300±200 


5.8±0.2 


6 


SMMJ1235+6215 


64.1 ±22 


0.63±0.06 


<4.1 


30.4±5.8 


37.4±6.3 


27.0±6.8 


8.3±2.5 


6.9±1.0 


1200±170 


0.7±0.4 


7 


SMM J 1237+6203 


(1.02±0.09)t 


(1.61±0.7)t 


<6.0 


<14 


<24 


<24 


5.3±1.7 


2.4±0.2 


430±50 


5.33±0.05 


8 


SMM J 1636+4057 


66.8±6.8 


0.95±0.1 


<9.1 


44.4±5.8 


36.6±6.3 


29.0±6.8 


8.2±1.7 


7.9±1.5 


1400±2()0 


1.9±0.6 



Table 2. Notes: 

Colum n (1): Source ID. Column (2): Source name. Column (3): 8pm flux densities and t h eir uncertainties take n from: Hainline et al. lfe009t);ICardamone et al] 
or the Spitzer archive. Column (4): 24pm flux densities from: llvison et al] f2007); Valiante et al. ( 2007); Menendez-Delmes tre et al l j2009|): Dickinson 
et al. i n prep; or the Spitzer archive. Uncertainties are assumed to be 10% of the measured values. Column (5): 70pm flux densities taken from lHainline et al] 
(2009). Columns (6)-(8): The 350pm, 450pm and 500pm flux densities. This HerMES data was accessed through the HeDaM database (http://hedam.oamp.fr) 
operated by CeSAM and hosted by the Laboratoire d'Ast rophysique d e Mars e ille. The uncerta i nties quoted are the confusion noise for each band. Co lumn (9): 
850 pm flux densities and their uncertainties taken from: ICowie et alj j2002h : IChapman et al] <2004 : IChapman et alj i2005h : IClements et al] (2008). Column 
(10): Derived IR luminosities (8-1000pm) from the star formatio n component of th e SED fitting procedures described in [[5] Column (1 1): Star formation rates 
derived using the quoted Lir.sf values and the relationship from Kennicu^ il998h . Column (12): The estimated AGN bolometric luminosity calculated using 
the methods outlined in j3] JThe value quoted is calculated using the radio-infrared correlation (assuming a q = 2.1). This is an upper limit as there will be 
an unknown AGN contribution to the radio flux. **This AGN bolometric luminosity is derived from the [O III] luminosity (see [[5). * For SMM J0943+4700 
the flux densities have been corrected for an amplification factor of 1.2. We note that the quoted flux densities will cover both radio counterparts to this source 
(see Appendix A5). tThese are the 12pm and 22pm flux densities taken from the Wide-field Infrared Survey Explorer archive. 

Galaxy-integrated [O iii]^5007 Emission-line properties 



Component A. Component B 



ID 


Source 


Z[oiii],a 


FWHM A 


FIuxa 


Z[om],B 


FWHM B 


FIuxb 


Av 


Total [O III] Luminosi 






(kms~') 


(10~ 16 ergs-' cm" 2 ) 




(kms^ 1 ) 


(Kr 16 ergs"' cm" 2 ) 


(kms" 1 ) 


(10 42 ergs"') 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1 


SMM J0217-0503 


2.02 13 [6] 


540±170 


3.1±0.7 










9.3 


2 


SMMJ0302+0006 


1.4049[10] 


900±300 


1.3±0.4 










1.6 


3 


RGJ0302+0010 


2.23861[12] 


200±40 


0.87±0.17 


2.2419[7] 


1100±120 


2.7±0.4 


±300±60 


13.9 


4 


RGJ0332-2732 


2.3145[5] 


480±140 


2.3±0.5 










9.9 


5 


SMM J0943+4700* 


3.3506[4] 


430±110 


1.1±0.5 


3.3503[8] 


1300±200 


2.6±0.4 


±20±60 


38.4 


6 


SMMJ1235+6215 


2.199[2] 


1500±500 


1.2±0.4 










4.5 


7 


SMM J 1237+6203 


2.07502[3] 


234±10 


4.1 ±0.2 


2.07 14[4] 


1000±100 


2.9±0.3 


-350±40 


22.5 


8 


SMM J 1636+4057 


2.38500[12] 


240±30 


1.3±0.3 


2.3817[18] 


900±200 


1.0±0.4 


-290±160 


16.8 












2.4003[10] 


1100±180 


1.4±0.2 


+1350±90 





Table 3. Notes: 

The properties of the galaxy-integrated [O III] emission-line profiles shown in Fig. [5] Numbers in square brackets give the uncertainty in the last decimal place. 
Column (1): Source ID. Column (2): Source name. Columns (3)-(5): Properties of the narrowest Gaussian components (Component A): redshift, FWHM and 
flux respectively. Columns (6)-(8): Same as the previous three columns but for the broader Gaussian components (Component B). Column (9): The velocity 
offset between the two Gaussian components. Column (10): The total [O III] luminosity. The bottom line of the table gives the properties of the third Gaussian 
component which is fit to the profile of SMM J1636+4057. Flux and luminosity values are not extinction corrected. The quoted uncertainties are from the 
parameters of the emission-line fits. The true uncertainties on the fluxes are a factor or «2-3 higher, due to uncertainties in the absolute flux calibration. *The 
flux values for SMM J0943+4700 have been corrected for an amplification factor of 1 .2. 



Spatially Integrated Emission-Line Properties 



to 
O 



[O III] Component A. [O III] Component B Hp 





ID 


Source 


Region 


Z[om],A 


FWHM A 


FIuxa 


Z[onrj,B 


FWHM B 


FIuxb 


Av 


'hp 


FWHM Hp 


Flux H p 


log([0 m]/HP) 


LinearSxtent 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(T5) 


1 


SMM 10217-0503 


BR 


2.0229[11] 


900±400 


0.42±0.15 








-460±60t 






<0.31 


>0.13 








NR 


2.0209[2] 


290±80 


1.26±0.19 


. 


. 


_ 




. 


_ 


<0.34 


>0.57 


n st 4 


3 


RGJ0302+0010 


BR 


2.23873[7] 


180±20 


0.53±0.06 


2.2410[3] 


1180±80 


2.59±0.16 


+210±30 


2.2441 [11] 


780±170 


0.34±0.10 


0.97±0.13 


8 6 






NR 


2.23864[12] 


200±30 


0.17±0.02 














<0.03 


>0.73 


6*£ 


4 


RGJ0332-2732 


BR 


2.3 143 [4] 


560±130 


0.81±0.16 








-100±20t 






<0.20 


>0.61 


10&4 






NR 


2.3147[4] 


510±130 


0.62±0.12 














<0.19 


>0.52 


ir±4 


5 


SMM 10943+4700* 


BR 


3.35 11 [4] 


340±70 


0.44±0.13 


3.3503[5] 


1400±100 


2.41±0.17 


-50±40 


3.3500[7] 


430±120 


0.30±0.09 


0.97±0.12 


15±3 






NR 


3.3493[2] 


400±40 


0.48±0.04 














<0.07 


>0.71 


8±3 


8 


SMM 11636+4057 


BR 


2.3834[8] 


1100±170 


0.77±0.12 


2.4023[7] 


1220±160 


1.07±0.14 


+ 1670±100tt 






<0.11 


>1.2 


7±3 






NR 


2.38481 [7] 


257±14 


1.25±0.06 














<0.11 


>1.1 


20±3 



Table 4. Notes : 

The properties of the [O III] and Hp emission-line profiles extracted from the sub-regions shown in Fig. [5] Column (1): Source ID. Column (2): Source name. Column (3): The spatial region from which the 
spectra are extracted, either the broad region (BR) or narrow region (NR); see i]4.3l Columns (4)-(6): Properties of the narrowest [O III] Gaussian components (Component A): redshift, FWHM (kms -1 ) and 
flux (10~ 16 erg s _1 cm -2 ) respectively. Columns (7)-(9): Same as the previous three columns but for the broader [O III] Gaussian components (Component B). Column (10): The velocity offset (kms~') between 
the two Gaussian components. Columns (11 ) — ( 13): Properties of the Hp emission lines: redshift, FWHM (tans -1 ) and flux (10~ 16 ergs cm ) respectively. Column (14): The logarithm of the ratio of the total 
[O III] and Hp fluxes. Column (15): The approximate linear-extent, along the most extended axis, of the defined regions (kpc). We take the size of seeing disks as the uncertainty on these measurements. All flux 
limits are 3o~ upper limits, which are calculated using the noise of the line-free continuum and the FWHM of the [O III] emission line from the same region. Flux values are not corrected for dust extinction. The 
quoted uncertainties are from the parameters of the emission-line fits. The true uncertainties on the fluxes will be a factor or s»2— 3 higher, due to uncertainties in the absolute flux calibration. *The flux values for 
SMM J0943+4700 have been corrected for an amplification factor of 1.2. t These are the velocity offsets to the Ha emission in this region, f t F° r analysis purposes we define the velocity offset as half of this value, 
see i]4.4| and Appendix A8 for details. 
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APPENDIX A: NOTES ON INDIVIDUAL SOURCES 

Here we provide background notes on each source from our sam- 
ple and provide some more details of our analysis on individual 
sources that is not required for the overall discussion given in the 
main text. The 1.4GHz flux densities are provided in Table[T]and 
the infrared and sub-mm flux densities are provided in Table [2] All 
of the parameters for the profile fits to the galaxy-integrated spectra 
(see Fig.[3} and spectra extracted from different regions (see Fig. [5} 
are provided in Tableland Table ^respectively. 



Al SMM J0217-0503 

This source was firs t ide ntified as a sub - mm a nd radio source by 
ICoppin et al] (2006) and ISimpson et al.l d2006l) respectively. IFU 
emission around the Ha emission line reveal at least two merg- 
ing components ( Alaghband -Zadeh et al ]|2012l ) over fa20-25 kpc. 
Our IFU observations of the [O III] emission line also reveal at 
least two systems separated by ~3" (~25kpc; see Fig. [4j. The 
northern system has broad Ha emission in the central regions 
(FWHM~940kms~ 1 ; Alaghband-Zadeh et al. priv. comm.) and 
a high [N IIl]/Ha ~ 1.1-1. 3 ratio (Alaghband-Zadeh priv. comm). 
Along with our measured ratio [O III]/HP >1.3 this suggests that 
this s ource is most likely hosting AGN activity (e.g. iKewlev et al.l 
2006). In Fig.[5]we show that the [O III] kinematics in this northern 
component also displays the signatures of an outflow; i.e. a broad 
emission line (FWHM=900±400kms- 1 ) which has a high veloc- 
ity offset from the Ha emission (Av = -460 ± 60 km s" 1 ; Fig. EJ. 
Due to the complexity of blending between the [N II] doublet and 
Ha line, it is not possible to confirm if this outflowing component is 
also present in the Ha emission-line profile. The broad [O III] emis- 
sion in this northern component is spatially unresolved (^ 4.0 kpc) 
and has a modest signal-to-noise ratio such that we do not have 
strong evidence for a galaxy-wide outflow in this source. 

The southern system has Ha emission extended over ~ 10 kpc 
with kinematics which may be described by a rotating disk; how- 
ever, the double peaked emission line and a large gradient in the 
[N Il]/Ha ratio from east to west could als o indicate that this south- 
ern sy stem has undergone a recent merger (Alaghband-Zad eh et al.l 
l2012h . Our observations of the [O III] emission line displays a red- 
shift, line-width and position that is only consistent with the west- 
ern ~5 kpc of the Ha emission, with no [O III] detected in the east- 
ern side. The [O III] emission is located where the [N Il]/Ha ratio 
is lowest suggesting a different chemical co mposition across th e 
galaxy, as opposed to excitation effects (e.g., lKewlev et al1l2006h . 
These observations are consistent with the scenario that this south- 
ern system may have undergone a recent merger. 
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A2 SMM J0302+0006 

This sub-mm source was first detected by IWebb et all d2003l) . 
Optical and near-IR spectra show a high [N Il]/Hoc =1.38±0.07 
line ratio, and Spitzer mid-infrared spectroscopic observa- 
tions show an exc ess in the mid-infrared continuum suggest- 
ing AGN activity JSwinbank et all |2004 IChapman etal] 120051 : 
Menendez- Delmestre et al.ll2009h 7 HST imaging and Hot IFU data 
show this source cont ains two components separated by ~ 1 1 kpc 
dSwinbank et al J2006I : Menendez-Delmestre et al. 2012). 

The galaxy-integrated [O III] emission we observe is broad 
(FWHM=900±300kms~ 1 ; Fig. [3) but is spatially unresolved « 
3.6 kpc). However, due to the low signal-to-noise ratio of the in- 
tegrated spectrum we cannot rule out the existence of low surface 
brightness extended emission. No Hp emission was detected in the 
galaxy-integrated spectrum, with a 3a upper limit on the flux of 
F H p < 8 x 10~ 17 erg s~ 1 cm~ 2 (assuming the line- width is the same 
as that measured for [O III]). 



appear to be tracing similar kin e matic s to the Ha emis- 
sion line ( Alaghban cUZadeh et al.l l2012h with similar line 
widths (FWHMss200-700kms _1 ) and velocity gradients 
(Av <200kms" 1 ). The [O III] emission is more extended to 
the north than the Ha emission and, along with t he observation that 
the [N Il]/Ha ratio is highest in the north (Alaghband-Zade h et al.l 
|2012|) . suggests the presence of an AGN ionisation cone. 

Unlike the majority of our sample with spatially resolved data, 
we are unable to fit multiple components to the emission line pro- 
file and are unable to de-couple the signatures of an outflow from 
the galaxy dynamics. However, we note that the emission lines 
in the brightest regions display very large widths (FWHM up to 
R*700kms , Fig. [6jl, which are broader than the narrow lines 
we associate with galaxy dynamics in the other sources (~100- 
500kms _1 ; Fig. |7](. Although not as well defined as in the other 
sources, we speculate that an outflow could be responsible for these 
large line widths. 



A3 RG J0302+0010 

This source was first identified by Chapman et al. (2004). This 
source has rest-frame UV spectroscopy that reveals strong [C IV] 
emission, which along with the the high [N H]/Hq=1.13±0.4 line 
ratio implies AGN activity dChapman et al. 2004; Swinban k et al.l 
120041) . Observations with long-slit spectroscopy show that the 
O III] emission is broad and has an asymmetric line profile 
Takata et al1l2006h . 

The galaxy-integrated spectrum around the 
[O III]XX4959,5007 emission-line doublet reveals both a narrow 
and redshifted broad component (Fig. [3}- The broad component is 
located in the central ~8 kpc, where the line ratio [O IIl]/Hp=9±3 
implies bright AGN activity (Fig. [5}- We confirmed the broad 
emission line was spatially extended by integrating spectra in 
two 0.6" x 0.6" bins, one in the north of the emission-line region 
and one in the south. We fit these spectra following the methods 
outlined in ij3]and found that the broad emission line has a velocity 
shear of 125±30kms _1 which is broadly consistent with the 
velocity gradient observed in Fig. [6] 

The broad [O III] emission line is most likely due to an out- 
flow with the near-side of the outflow being obscured by dust (see 
£]5.2| for a full discussion). Blueshifted br oad components are m ore 
commonly associated with outflows (e.g. jHeckman et alJI 198 lb as 
opposed to the redshifted broad component we observe here. How- 
ever, redshifted broad components are observed in a fraction of lo- 
cal AGN and can be explained by the relative orientations of the 
AGN and obscuring m aterial in host galaxy (e.g jBarth et al .120081 : 
ICrenshaw et al.ll2010l) . In particular IFU observations of the local 
interacting Seyfert galaxy LEDA 135736, identify a broad, red- 
shifted outflow associated with AGN activity, which could be a low 
redshift analogue of this source jGerssen et ai]|2009l) . 



A5 SMM J0943+4700 



This sub-mm source was first identified by ICowie et al.l (2002) 
and is modestly lensed (amplification of 1.2). lLedlow et al.l (2002) 
identified two radio counterparts separated by ~30kpc in pro- 
jection and labelled them as H6 and H7. The latter shows nar- 
rower and fainter [O III] emission than its companion but dom- 
inates the sub-mm and CO emission, indicatin g that the bulk of 
star formation is oc c urring in this component (jj akata et al. 2006; 
IValiante et ail 120071 ; Enge! et alj |2010| ; iRiechers et al.ll201ll) . An 
excess in the mid-infrared continuum observed in Spitzer mid- 
infrare d spectroscopy indica te the presence of AGN activity in this 
source dValiante et al]|2007h although both H6 and H7 are covered 
by this spectra. H6, the source observed in our observations, has 
rest-frame UV emission-line wi dths, and line ratio s indicative of 
a narrow-line Seyfert 1 galaxy dLedlow et alj|2002l) . Additionally 
long-slit near-IR spectroscopy of H 6 reveals [O III] em ission that 
is broad and extended over as3"-4" ( Takata et al. 2006), indicative 
of an outflow. 

Using our IFU observations we find that the [O III] emission 
from this source is very extended (~20kpc) with a complex mor- 
phology and kinematic structure (Fig. [4}. The north-western region 
shows modest [O III] line- widths (FWHM=200-300kms~') con- 
sistent with being due to galaxy dynamics and merger remnants 
(Fig. [7} see § 15. II for a discussion). In addition there is extremely 
broad [O III] emission (FWHM= 1000- 1400 km s _1 ) extended over 
~15 kpc (Fig.[5]and Fig.[6}, which we attribute to an energetic out- 
flow. Despite observing modest velocity offsets between the narrow 
and broad [O III] emission lines (Av ^150kms _1 ; Fig. [6]», there is 
clearly extremely turbulent ionised gas over a considerable extent. 
The zero velocity offset between the broad and narrow [O III] com- 
ponents in the spatially-integrated spectra (Fig.|3]and Fig.[5j( could 
be explained if we are observing an outflow orientated in the plane 
of the sky. 



A4 RG J0332-2732 

IFU observations of this source reveals that Ha is extended over 
~ 10 kpc with velocity offsets of Av < 200 km s~ 1 across the galaxy. 
The brightest region (Fig. [5] Table g} has a high [O III]/H(3 >4.0 
line ratio, which along with t he high [N Il]/Ha=0.8±0.1 
line ratio (A laghband-Zad eh et all 1201 2|) suggests AGN activity 
(e.g jKewlev et alj200d) .~ 

The [O III] emission-line kinematics (Fig. |4j Fig. [6} 



A6 SMM J1235+6215 

This source was id entified as a sub-mm bright galaxy by 
Chapman et al. (2005 ). Deep Chandra images shows that it hosts a 
heavily obscured X-ray luminous AGN (Lo.5_8keV = 10 44 ' erg s _1 ; 
Nu=10 cm " 2 ; lAlexander etal] [2005). An excess in the mid- 
infrared continuum observed in Spitzer mid-infrared spectroscopy 
further suggests the presence of AGN activity in this source 
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dMenendez-Delmestre et aD |2009). IFU observations reveal a 
bright compact (~3kpc) source of broad Ha emission (FWHM 
> 1000 km s -1 ) attributed to AGN activity and spatially off- 
set (fa0.5") from the narrower Ha emission (FWHM 500- 
800 km s ), which is likely to be dominated by star-formation ac- 
tivity (Menendez-Delmestre et al. 2012). Long-slit near-IR spec- 
troscopy su^g£stjhatthe_[0 III] emission is extended over fa 1" 
(fa8 kpc; lTakata et alj2006h . 

The galaxy-integrated spectrum around the 
[O III]A,A4959,5007 emission-line doublet reveals an extremely 
broad emission line (FWH M= 1 500±500 km s~ 1 ) which may 
suggest an outflow (Fig. [3). iTakata et alj d2006t) find that this 
sources has extended [O III] emission over fal". Although we find 
tentative evidence that the [O III] emission is marginally extended 
(fa0.8"; 7 kpc) in our data we are unable to constrain the properties 
of the extended emission due to insufficient signal-to-noise ratio of 
the data to fit the emission-line profiles. 

We also detect Hp emission, at a redshift of 
z H p =2.203 1 ±0.0006, with a width of FWHM H |3=500±200 km s" 1 
and flux F H |3 = 0.9±0.2x 10~ 16 erg s _1 cm~ 2 , extended over fal" 
(~8kpc). The extremely broad [O III] emission we observe is 
offset from the Hp emission Av fa -400 ± 200 km s" 1 . Although 
there is weak evidence for an extended outflow in this source, 
deeper observations are required to confirm this. 



A7 SMM J1237+6203 



This source was ide ntified as a sub-mm bright galaxy by 
IChapman etail d2005h . The bright optical counterpart (R = 
l2005h. 



20.2; 
10 44.3 



Chapman et al. 



X-ray luminosity (Zo.5-8keV = 
erg s~ l ; lAlexander et al.1 [2005b and optical-near-IR spec- 
troscopy that reveal broad emis sion lines (Lyq; N V; C I V; Hy-Hq; 
FWHM fa 2100-2700km s~ 1 ; IChapmanetalJl2005l :l Takat a et al.1 
l2006l : ICoppin~ et alj|2008l) . result in this source being classified as 
a broad-line quasar, w ith a virial BH mass of log(MBH)=8.2MQ 
( lAlexander et alj|2008h . 

IFU observations around redshifted [O III] of 
SMM J1237+6203 were first presented in lAlexander et al.1 
where a complete discussion of the source is pr o vided . 
We note the orientation of this source in lAlexander etail j20 id) 
is flipped in the east-west direction and the correct orientation is 
given in this work (Fig. [4). In addition we provide updated [O III] 
flux and luminosity measurements. 

The galaxy-integrated spectrum around the [O III] emission 
line (Fig. [3j shows a bright narrow component with a prominent 
blue-wing. This type of profile is most comm only interpreted as be- 
ing the result of an AG N-drfven outflow (e.g., Heckma n et al.ll98ll ; 
INesvadba etai]|2008h . The velocity and FWHM profiles (Fig. [6} 
show that the narrow emission has a small velocity gradient of 
Av fa 200 km s~' across the ga laxy and is likely to b e tracing the 
host galaxy dynamics (also see Alexander et al.1 2010l) . In contrast, 
the broad [O III] component (FWHM fa 1000 km s~ 1 ) is offset from 
the narrow component by fa — 350kms _1 , over ~5kpc, providing 
strong evidence for an outflow in this region. Although the broad 
component is only formally fit over ~5 kpc, low s urface brightness 
broad emission exists up to ~8 kpc in total extent (Alexan der et al.l 

12013). 

We also identify HP emission in the galaxy-integrated 
spectrum, with a redshift of z=2.0754[2] and flux 
F H(3 =l.l±0.4xl0" 15 ergs- 1 cm- 2 . The width of the Hp 
emission line (FWHM=2030±70kms~') indicates the presence 



of a broad-line region, which is consistent with that found from 
Ha observations of this source dCoppin et al.ll2008h . 



A8 SMM J1636+4057 



This sub-mm so u rce w as first identified by IScott et al l J2002l) 
and llvison et alj d2002l) . Hubble Space Telescope (HST) ACS 
and NICMOS observations reveal a complex galactic morphol- 
ogy with at least three mergin g or interacting com ponents, cov- 
ering ~8-15 kpc in projection dSwinbank et alj2005t see Fig. lAlt . 
An excess in the mid-infrared continuum observed in the Spitzer 
mid-infrared spectros copy indicate the presence of A GN activ- 
ity in this source ( Men endez-Delmestre et al] [2009). Long-slit 
and IFU near-IR spectroscopy reveal a broad-line AGN (broad 
spatially unresolved Ha and HP emission with FWHMfa2000- 
3000 km s" 1 ) and narrow Ha emission traces star-for mation along 
the UV bright arc shown in the HST image (Fig. IA11 ISmail et all 
l2003l : ISwinbank et al.120051 ; Menendez-Delmestre et al. 2012). The 
broad-line AGN is found at Av fa 700kms _1 from the narrow 
Ha emission, providing further evidence that this system is un- 
dergoing a merger. The presence of a molecular gas cloud in the 
southern regions of the system i s shown through h i gh and low 
excitation CO emission-line gas dTacconi et al]|2008l ; Enge l et al.1 
l201Ctllvison et alj|20l7l) . Previous long-slit observations show that 
the [O III] emission is kinematically complex and spatially ex- 
tended (fa 1-2") with broad compon ents (FWHM < 2000 km s _1 ; 
ISmail et alj2003l ; lTakata et al.l2 006). providing initial evidence for 
an energetic outflow. 

Our galaxy-integrated spectrum around the 
[O m]M4959,5007 emission-line doublet (Fig. O reveals a 
prominent narrow component and two broad components sepa- 
rated by Av = 1700 ± 100 km s" 1 . We are confident that the broad 
redshifted emission line is also from [O III] as lSmail et al.1 (2003) 
also identify [Ne V] and N V emission lines at approximately this 
redshift. 

In Fig. [4] we show that the [O III] emission-line morphology, 
velocity field and FWHM map from our IFU observations. Narrow 
[O III] emission (100 < FWHM < 500 kms" 1 ) is extended over 
fal" (fa8 kpc) to the north of the peak in surface brightness, associ- 
ated with faint infrared emission in the HST image, and fainter nar- 
row emission is observed fa l" to the south; th is is consistent with 
the long-slit observations o f lSmailetal.l l l2003l) (also see Fig. lAU . 
The redshift of this narrow [O III] emission is consistent with the 
molecular gas and Ha arc and therefore is likely to be tracing 
merger remnants and galaxy dynamics (e.g., Fig.|6). We also detect 
faint HP emission that is spatially coincident with the star-forming 
arc (Fig.lAlY 

The two broad [O III] emission-line components 
(FWHM~1200km s _1 ) are spatially coincident with the broad Ha 
emission line (Fig. [5] and Fig. IAU and are moderately spatially 
extended (~0.8"; ~7kpc). We verified that the broad emission 
lines were intrinsically spatially extended by extracting spectra 
from two 0.6" x 0.6" bins; one in the north and one in the south 
of the observed broad [O III] region. We found that that the 
ratio of the two broad components changed from north to south. 
The redshifted broad component dominates in the north and the 
blueshifted broad component dominates in the south, which is also 
demonstrated in the velocity map of this source (Fig.|4). The broad 
[O III] components are offs et from the broad Ha emission line by 
Av fa ±850km s _1 (Fig. |U ISwinbank et alj|2005l) , indicating that 
we ar e observing both sides of an outflow (e.g., INesvadba et al.1 
2008). To verify that the broad Ha is due to a true broad-line 
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Figure Al. Left: True colour / 8 i4#i60 image of SMM J 1636+4057 from the H ST ACS and NICMOS imaging (Swinbank et al. 2005). The contours denote 
the intensity of the velocity integrated CO(7-6) millimetre emission iTacconi et alj2008t) . The circle denotes the approximate seeing disk of our observations. 
Middle: HST NICMOS H\(® grey-scale image with contours overlaid of narrow band (wavelength collapsed) images from our IFU observations. The solid 
contours denote the extent/intensity of the broad [O III] emission, which is associated with the red continuum seen in the colour image. The dotted contours 
denote the extent/intensity of the narrow [O III] for which the majority is associated with faint emission observed to the north of the red continuum peak. The 
dashed contours denote the Hp emission which is spatially coincident with the blue arc in the HST image. Right: The same image as in the central panel but with 
the contours from the IFU+adaptive optics Ha observations of Menendez-Delmestre et al. (2012) overlaid. There is extremely broad (FWHMai2700 km s ) 
unresolved Ha emission, attributed to an AGN broad line region (BLR) and narrower spatially extended emission along the blue arc, which is attributed to star 
formation. The field-of-view of these Ha observations did not cover the region of narrow [O III] emission seen in the north of the middle panel. North is up 
and east is left. 



region, as opposed to d ynamics in the host ga laxy, we attempted 
to re-fit the spectrum of ISwinbank et ay (120051) with a narrow-line 
Gaussian component plus two additional Gaussian components 
with the same redshifts and line-widths of the br oad [O III] emis- 
sion lines (e.g., similar to the methods used by ICano-Dfaz et al.l 
120121 and Mullaney et al. in prep). We found that this did not 
adequately describe the Ha emission-line profile and an additional 
broad component is still required (with FWHM~2700 km s -1 and 
z = 2.393 ±0.005; consistent with Swinb ank et alj2005l) providing 
strong evidence of a true AGN broad-line region. 

The outflow we observe is spatially coincident with the AGN 
in this source, and no outflow is observed in the region of unob- 
scured star-formation (the UV bright arc). We conclude that the 
broad [O III] emission lines we see are due a bi-polar AGN-driven 
outflow at v ~ ±850kms _I . 

A9 SMM J2217+0010 

This source w as first identified by ISmail et al l d2004 and 
IChapman et ail ( 120051) . An excess in the mid-infrared 
continuum observed in the Spitzer mid-infrared spec- 
troscopy sugge sts the presence o f AGN activity in this 
source (Menendez-Delmestre et al. 120091) . The results 
of iTakata et al.l J2006I) suggest bright [O III] emission 
(F [om] =6.0 x lO^^ergs- 1 cirT 2 ) and a high [O III]/HP «10 
ratio further indicating AGN activity. However in our observations 
[O III]AX4959, 5007 emission was undetected. By spatially inte- 
grating over the extent of a seeing disk and assuming an emission 
line with a FWHM of 500 km s , we infer a 3o upper limit on 
the [O III] flux (luminosity) of F[ 0m j <4.6 x 10~ 17 ergs -1 cm~ 2 
(Ljom] < 2.6 x lO^ergs -1 ) . This upper limit is considerably 
lower than the value quoted in lTakataetal.N2006l) which may arise 



due to the presence a strong sky-line at the observed wavelength of 
the [O III]A,5007 emission line in this source, potentially resultin g 
in an excess measurement of the flux quoted in Takata rtakl d2006r) . 

A10 SMM J2217+0017 

This source was i dentified as a sub-mm bright galaxy by 
iTamura et all J2009I) - We did not detected [O III]M4959,5007 
emissi on. This source was detected in Ha emission in the same data 
cube flAlaghband-Zadeh et al.l2012l ), therefore the non-detection is 
due to an [O III] flux that is below our detection threshold. Assum- 
ing an emissi on line of a FWHM of 350 km s _1 (to match the Ha 
observations; Alaghband-Zadeh et al. 1 120121) and integrating over 
the spatial extent of the Ha emission, we infer a 3c upper limit on 
the [O III] flux (luminosity) of Fiorrr] < 2.6 x 10~ 17 ergs~' ctrT 2 
(L|om] <1.0xl0 42 ergs- 1 ). 



